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The resurrection and the success of activities addressing all flight regimes depend critically upon the evaluation and 
understanding of the viscous flow effects. In the supersonic and hypersonic regimes, the boundary layer status affects 
both the friction drag and the thermal loading of the structure, two items of high consideration in any design. In the sub- 
and transonic regimes, sustaining laminar flow over a significant part of the wings, is a major challenge of current 
aerodynamic interest. ' 

In spite of the significant potential demonstrated by IR imaging systems for aerodynamic research, the method was not 
widely accepted and recognized for its full capability until relatively recently. Two problems delayed acceptance of the 
method in its earlier stages of development. The first was a lack of easily accessible capability to store and manipulate 
data. This drawback was naturally solved with the commercial introduction of video cassette recorders (VCR's) and 
personal computers. The second was the realization that the temperature information contained on thermograms results 
from a multi-faceted process encompassing aerodynamics, heat conduction, geometry, and radiation. 2 3 

The reduction of the experimental data, from temperatures to heat flux distributions via the thermal response of the 
substrate, can be simplified if the the physical reality can be mathematically modeled based on assumptions leading to 
analytical solutions. At one end are models made of thermal insulating materials, where it is plausible to assume that 
deep m the substrate there is a layer that virtually keeps its temperature constant throughout the test (if its duration is not 
too long). This case is known as the "semi-infinite slab" model, for which the heat conduction equation can be solved 
analytically for the heat flux as a function of the transitory surface temperature distribution. At the other end there are 
hollow, "thin skin" models, where it can be assumed that the heat flux is uniformly absorbed in-depth the material with 
the temperature changing uniformly across the model skin. This case can be approximated by the "thin skin" model of 
the thermal energy equation. Other types of model construction may fall in between these two classifications, as they are 
affected by the filler or the spar underneath the skin. It follows that understanding the direct and indirect factors 
influencing an experiment and its results, is critical for the determination of the data reduction technique. Therfore, 

designing aerodynamic experiments based on IR imaging and extracting data from the thermograms are very challenging 
processes. 66 

Utilization of the technique has yielded a wealth of experience and aerodynamic results, but the subject has not been 
adequately reviewed. As recently as 5 years ago, a paper on developments in flow visualization 3 mentioned only two 
references about aerodynamic research performed with IR imaging systems, although by that time the method was known 
for 21 years. This paper is a wide ranging review of the applications IR imaging has seen in aerodynamic research It 
should benefit m particular those contemplating the use of the method. After a brief description of IR imaging systems, 
a few overview papers will be mentioned describing and analyzing the IR imaging technique for aerodynamic research. 
Attention will turn then to super- and hypersonic, and in turn sub- and transonic applications, two areas where these 
systems made distinctive contributions. The survey will then proceed to some specific applications in flow visualization, 
propulsion, fluid mechanics and heat transfer. Before concluding, a quick look will be taken at current trends and 
expected developments in the application of IR imaging to aerodynamic research. 


2, IR Imaging Systems and Data Processing 

Most of the commercial IR imaging systems are similar in concept to those built 25 years ago. A single detector, usually 
cooled by liquid nitrogen, is exposed to the incoming radiation from the target on which the objective lens is focused. 
To obtain an image of the objective lens field of view (FOV), the camera scans the scene both horizontally and 
vertically, with rotating prisms or mirrors, producing the display lines for the video display. The area projected at any 
single time on the detector is subtended in the instantaneous field of view (IFOV). The output of the detector is 
transmitted to a black-and-white video display, where lighter shades are associated with higher temperatures. Usually, a 
single pixel displays the sensor output from one IFOV. To give some specific numbers, 4,5 commercial imagers are 
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targets and vice versa. The field-of-view of commercial objective lens ranges between 2.5 o g * P 

targets and vice ^versa ^ . ve the visua i quality of the display, two or four fields are usually 

r^teriaced°to produce a frame. Thus, the camera produces fields, but the human eye perceives fram ®*_ A “ 

usually fewer than 200 lines, each line having usually less than 200 pixels. The scanning rate may _vary A 

fields per second. Gauffre and Fontanella. 6 give particularly concise and uselul description of IR imaging y. 
introduces the reader to camera figures of merit, including sample calculations. 

_ . . tn vrR „ for analog recording and play-back and to personal computers equipped with 

evaluated on an input-versus-output basis. 9 In other cases, specific perfomumce £ 

Function (MTF) of the system may have to be evaluated, to account for mtsstng mformatton. and allow data 

restoration. 10 

™ sidltions, environmental comp.ltbility, and cost, Antimflective coating may improve the tmnsnutunce of some of 
these materials, but at the expense of the reflectance. 
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4, Super- and Hypersonic Studies 


The first documented use of IR imaging in aeronautical research was reported in 1967 by Thomann and Frisk 24 of the 
Aeronautical Research Institute of Sweden. In a wind tunnel experiment performed at Mach = 7, the temperature 
distribution on an elastomeric paraboloid was measured as it evolved in time. From the temperature rate of change the 
heat flux at the model surface was deduced using the "semi -infinite slab' solution 25 for the unsteady heat conduction 
equation. The short test run (seconds) combined with the low thermal diffusivity material justify use of that solution. 
This data reduction concept was borrowed from the fusible paint technique, 26 and it is still in use. The experiment 
showed the IR imaging to be as accurate as other competing techniques for heat flux determination, with the advantage 
that without requmng surface preparation before or in-between runs, it is quite expedient. More than a decade later, 
Balageas and Ory proposed a data reduction method that incorporated the finite thickness of the skin and the boundary 
conditions at the internal wall, expanding the applicability of the technique to thin-skin models. 


After this technique was shown to work, there were more attempts to evaluate and improve it, e.g., scanning the model 
against a water cooled plate to reduce the background radiation noise. 28 However, for a few years the method was far 
from being productive enough to be directly applicable to project designs. Compton 29 at the NASA Ames Research 
Center realized that the bottleneck of the technique was the data acquisition, storage, and processing. One should realize 
that the heat flux distribution is calculated from the temperature readings on a pixel-by-pixel basis that was generated at 
rates of approximately 88000 IFOVs per second. The solution was devised to record the data on an analog tape, to 
digitize and read it into a computer for processing, and thereafter to display and plot the results. This engineering 
concept set the pattern for similar systems to this very day. Automating the data processing system gave access to all the 
information produced during a test. In particular, it was understood that monitoring the history of the heating rate for 
relatively long periods of time, could indicate if the boundary-layer transition front moved on the model during the test. 

In l9 J?' Arnold Eng,neenn S Development Center (AEDC) embarked on a large-scale program to develop the 
capability of extensive heat transfer testing in the hypersonic regime with an IR imaging system. 30 * 1 The von Karman 
facility was designated for that purpose, and engineering modifications were made to allow hosting an IR imaging system 
for test series that extended over long periods of time. An automated data processing system was developed that accepted 
input from the IR camera, thermocouples on the model, and temperature and pressure probes providing the operational 
parameters of the wind tunnel. In parallel, materials used for model manufacture were screened for compatibility with 
the technique and the environment, and their thermal and radiative properties were documented. Accurate numerical 
values of these properties are critical for processing the output data from the IR imager, and deducing of the temperature 
distributions on the model. To assess the accuracy of the technique, calibration procedures were developed; the 
repeatability of the measurements was evaluated; and, a measurement error model was implemented. 32,33 The results 
were confirmed also by other experimental means and were found to agree well with theoretical predictions. On the 
negative side, it was found that step changes in the temperature distribution could not be resolved, the response being to 
smear* or blur the result over a few adjacent pixels. Moreover, the camera displayed a consistent measurement error 
when the target temperature gradient exceeded a certain value. These shortcomings were blamed on the relatively large 
size of the imager's IFOV, but this aspect was not pursued enough to a full understanding. While a smaller IFOV will 
allow a better capture of abrupt changes in the temperature distribution, the blur effect will always occur, because of the 
Optical Transfer Function (OTF) of the camera. 

The infrastructure and expertise developed at AEDC in the 1970's was used to measure convective heating rates on a 
0.04-scale model of the Space Shuttle orbiter under flow conditions prevailing during the re-entry phase (Fig. 1) The 
test flow conditions were typically: Mach = 8, temperature 1300 R, Reynolds number 0.5X10 6 to 3.5X10* per foot, 
angle of attack 30 to 45 degrees. The data obtained were incorporated in the design of the thermal protection system of 
the orbiter. ’ Besides project-oriented applications, general heating studies of elastomeric materials were carried 
out, the results giving insights to phenomena that were harder to access with other diagnostic tools. For example, in a 
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testing ptoenun nin at Mach- 10, the surface heating was observed to be oonunifonn to a greater extent than could have 
forLit Hot streaks that may be have been caused by wing-tip crossflow effects were detected and merged 
Sing rates were obsetved when bubbles caused by local failure of adhesive bonds formed m the substrate. Even 
though part of the data could not be reduced using the "semi-infinite slab” solution, the underlying assumptions emg 
no longer valid, some of these effects could still be quantified, and operational conclusions could be drawn. 

Tests under rarefied flow conditions at very high Mach numbers, say Mach=20, are characterized by a short duration of 
a very few seconds and heat fluxes that may get as low as 0.5 kW / m 2 . In these cases, solid models may no longer t* 
acceptable because the relatively long time required for the temperature pattern to become established ,s unavailable For 
this application thin skin models and a conforming date reduction model may be the only way to get useful results - 
£lagSs et al 1215 extended the use of thin skin models to low Mach numbers (down to Mach -2.0) through stimulated 
fheXmphy". According to this concept, the model in the wind tunnel ,s heated by a bank of lamps over a pre-set 
period offime, and the heat transfer coefficient is determined from the temperature response of the m< ^ el in ' eg ™ ^ ° “ 
a period of time. This approach is less susceptible to curvature effects of the surface and noise, the date reduction being 

based on integration and not on differentiation. 

Next to determination of heat fluxes, location of boundary layer transition to turbulence was always a major subject of 
interest to aerodynamicists. Turbulent boundary layers have higher skin fnction than the laminar ones. On the o 
hand the laminar boundary layer is more susceptible to separation, especially in the super- and hypersonic regimes 
where shock-wave interaction can cause strong adverse pressure gradients. With the increase in skin fnct.on, the 
transition to turbulence also induces h.gher heat fluxes, and under some speed and flight time combination, the structure 
may require special thermal protection. Therefore, the capability to detect boundary layer transit, on may have 
far-reaching consequences for the development and verification of designs. 

The thermal signature of trans.t.on on a model can be caused by one of two effects. The first is the near tenn thermal 
response of the model, occurring as long as it does not reach thermal equilibrium with the recovery temperature in the 
horary layer. Under these circumstances, heat transfer takes place between the flow and the model, the area exposed 
tothetmbulent regime changing its temperature faster. 37 The second effect ,s the long term response, occurring when 
the model and the surrounding flow come to quasi-steady thermal equilibrium In this case, the wall tempera ure en 
towards the adiabatic value, the latter being higher under the turbulent-versus-lanunar boundary layer reg,me. 

Peake at al 38 carried out a boundary layer transition detection test in a blow-down tunnel at Mach number 3.85 using a 
stainless sted flat plate equipped with a bakelite insert surface. The transition location appeared on the thermograms as a 
localized hot frontthat Tas attributed to the difference in the recovery temperature between the laminar and the turbulent 
regime- it was also confirmed by other experimental means. In a subsequent experiment an all stainless stee flat plate 
painted black was used as target, but the distinct transition pattern could not be observed anymore.^ The inadequacy of 
Sie IR imaging technique to detect transition on stainless steel models is caused by the relatively high thermal diffusivity 
of tie m3 that leve.s-out temperature differences along and in-depth through the model. In cases where the surface ,s 
polished, the situation is ftirther aggravated. The high reflectance and the low em.ttance of the sur a “ ^r e ase he 
signal-to-noise ratio on the thermograms to an extent that renders the method useless. It may be speculated that g 
Mach numbers the very intense heat transfer may offset the detriments of polished stamless models by raising ery 
rapidly the surface temperature to relatively high values. Collier et al. 3 tried this approach on a polished stamless steel 
cone at Mach= 14 but, unfortunately no reliable temperature readings could be obtained. 

Identifying the indication of transition on thermograms with a precise locus inside the evolving transition Process is 
critical in evaluating the merit of the method. In more concrete terms, the problem is that on thermograms the transition 
to turbulence appears to happen abruptly, while actually the process may develop over a significant ^ portion of the 
aerodynamic surface. Hall et al. 40 compared the diagnos.s of transition from an IR imaging system with that of h 
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on a flat plate at Mach numbers 1.5, 2.0, and 2.5. The results indicate that on thermograms the transition appears at 
various stages of its development. As expected, it always appears before the 50% interim ttency; this interim ttency value 
decreases as the Mach number increases. 

Besides transition, flow separation and reattachment are features of much interest, especially on configurations of 
longitudinal symmetry. As the angle of attack increases, the flow wilt separate along the sides of the models, developing 
into large-scale vortices that reattach on the leeward side. These features are detectable through the surface temperature 
variation, minimum at separation and maximum at reattachment, and are induced by the respective behavior of the skin 
friction and heat transfer. Bandettini and Peake 41 carried out a separation detection study on a 10 degree fiberglass cone 
at different angles of attack at Mach= 1.8, comparing the results with oil-flow visualization. Although the IR imaging 
and the oil flow visualization were performed at somewhat different Reynolds numbers (3x 10 6 for the former and 9x10* 
for the latter, based on the model length), the vortices separation and reattachment could be identified on the 
thermograms as areas of higher and lower temperatures. An investigation of a somewhat more complicated flow pattern 
was carried out by Arai and Sato 42 on a 12.84/7-degree, bent-nose biconic made of epoxy resin at various angles of 
attack at Mach = 7. This geometry is a leading candidate for the forebody of the Aeroassisted Orbital Transfer Vehicle. 
Arai and Sato results are very similar to those of Bandettini and Peake, 41 confirming that vortical flow reattachment 
can cause heating on the leeward side at positive angles of attack. 

The experiments reviewed so far address temperature signatures of flow fields on models of circular geometry or 
longitudinal symmetry. Irately, Henckels and Maurer published the results of an experiment addressing a Mach = 8.7 
flow along a comer, produced by the intersection of two perpendicular plates, at various pitch and yaw angles. This is 
highly three dimensional flow, involving interaction of two perpendicular shock waves originating at the leading edge of 
each plate, and the formation of an embedded vortex underneath the slip surface produced by the interacting shocks. The 
thermograms were used to deduce heating rates on the walls, with the stipulation that the influence of the flow scaling on 
the results could not be inferred. 


Space Shuttle Flight Experiments 


The Space Shuttle program deserves a special place in the history of IR imaging in aerodynamic engineering and 
research. As it was previously mentioned, the method was used during the engineering phase of the program to 
determine the atmospheric re-entry heating rates on its forebody. After the orbiters became operational, they were 
designated for two daring IR imaging flight experiments aimed at mapping their surface heating during the actual 
re-entry phase. 


In the InfraRed Imaging of the Shuttle (IRIS) experiment, 44-46 the windward side of the Space Shuttle re-entering the 
atmosphere was observed from below through an astronomical IR telescope mounted on a 'chasing' C-141 aircraft, 
known as the Kuiper Airborne Observatory. In spite of useful and promising results regarding the evolvment of the 
actual heating rates and the progression of the boundary-layer transition front, 47 this experiment was discontinued 
because of its cost, the complex coordination work, and restricted availability of the airborne observatory. 

In the on-going Shuttle Infrared Leeside Temperature Sensing (SILTS) experiment, an IR imaging camera installed atop 
the vertical stabilizer of the Space Shuttle Columbia (Fig. 2), is scanning the leeside of that vehicle during the 
atmospheric re-entry phase of its flight. 48 ’ 49 Zones of intense heating were observed along the leading edge and the 
upstream part of the wing, the inboard/outboard elevon gap, and the orbital maneuvering system pod (Fig. 3), atop 
which the vertical stabilizer is mounted. 50 ’ 51 In general, the heating pattern indicates a highly vortical flow, converting 
high temperature compressed air from the windward side of the orbiter to its leeside. This behaviour is typical to 
delta-like configurations at high angles of attack experiencing leading edge flow separation. The real pay-off from this 
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experiment will materialize when the heating patterns observed in flight will be reproduced w numerical simulations of 
the flow. Thereafter, the in-depth thermal response of the wing substrate could be analyzed in detail, promoting the 
understanding necessary for future thermal protection designs. 


5. Sub- and Transonic Studies 

The possibility of using IR imaging for sub- and transonic aerodynamic research was raised only in the early eighties. 7 
Since then, the application of this technique got to the point where some industrial wind-tunnels offer IR imaging 
diagnosis as part of their standard data acquisition process.' 2 The interest in thermography for these flow regimes is 
mainly for boundary layer transition research. The method is particularly attractive because it produces global views of 
the configurations of interest, mainly airfoils and wings, where the behaviour of the boundary layer can be deduced 
through a quick visual inspection. 

The preoccupation with laminar flow for drag reduction reappears periodically in the aerodynamic research. In the past, 
it was used mainly as a means for range extension. More recently, the global energy crisis has caused the subject to 
resurface. New designs of transport airplanes are required to reduce their fuel consumption and their resulting direct 
operating costs. One avenue that can be used to achieve this goal is to design the wings to sustain laminar boundary layer 
to the largest possible extent. In the current laminar flow research, the IR imaging technique is an integral part of the 
effort to make this aerodynamic technology available. The effort encompasses wind tunnel technology, wind tunnel 
and flight testing, 54 ' 56 and aerodynamic design philosophies. 

In 1983, Bouchardy, Durand and Gauffre 7 from ONERA Chatillon reported the first transition detection study in a 
thermally stable wind tunnel. They addressed the effect of positive step change in the adiabatic wall temperature induced 
by the difference in the recovery temperature of the turbulent- versus-laminar regime. Using models fabricated of 
thermally insulating materials, current IR imaging systems are sensitive enough to detect transition down to Mach =0.5, 
at ambient conditions. For lower air velocities, down to Mach =0.1, the transition detection is still feasible through 
image processing enhancement of the thermograms. Among the thermally passive approaches to transition detection, 
this is the only aerodynamic effect that can be used in open circuit or thermally stable wind tunnels, where the stagnation 
temperature is considered constant. 60 The first systematic investigation of transition detection in wind tunnel and flight 
testing was publisheded by Quasi from DFVLR Braunschweig He pointed out that the increased heat transfer 
coefficient of the turbulent-versus- laminar regime is a remarkably useful attribute that makes the transition detectable on 
raw thermograms even at very low speeds, if the temperature of the model is different from the temperature of the flow. 
When exposed to this effect, the area underneath the transitional and turbulent boundary layer will change its temperature 
faster than its counterpart under the laminar regime, thus enhancing the temperature contrast between the two areas on 
the thermogram. In fact, the heat transfer effect is dominating over the recovery temperature influence, and today it is 
used in the majority of transition detection experiments. It is convenient in use and interpretation, especially implying 
the Reynolds analogy between convective heat transfer and skin friction coefficients. This effect occurs in closed circuit 
wind tunnels, where the air heats naturally during their operation. 62 It can also be observed in wind tunnels equipped 
with a cooling system. 61 Another option is to use electrical 19>20>63 or laser 64 heating of the model surface, or to blow 
hot air through a hollow model, 17 ’ 1 * and to observe the cooling effect of the flow. In flight testing, this effect occurs 
naturally when an airplane is flying an ascent or descent path through the atmosphere. Although active heating of the 
substrate is very attractive and simple to implement, the surface overheat value should be kept at the absolute necessary 
minimum to prevent a premature triggering of transition. 

When the significance and potential of IR imaging became evident for boundary layer research, the effort turned towards 
a full evaluation of the diagnostic capabilities of the method. In an effort to make the method more quantitative, it was 
shown that the active heating method can identify changes in the boundary layer regime by tracking the behavior of the 
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experimentally deduced heat transfer coefficient. Provided laminar boundary layer is established downstream of the 
leading edge, the maximum heat transfer coefficient is indicative of transition to turbulence, and its following minimum 
occurs at separation. 17 ' 20,63 ‘ 65 If the passive method is applied, the transition is identified visually at the location of the 
step change in the surface temperature. 7 ’ 12 ’ 55,61 ’ 62,66 ’ 67 At higher incidence, the onset of turbulent separation is 
identified by the emergence of a second zone of increased heat transfer caused by local vorticity shedding. 

The next step was to apply the IR imaging method to boundary layer transition research, that eventually would lead to 
"laminar" airplane designs. For this application, one must be able to identify on the thermograms the various modes of 
boundary layer transition, since the development of each mode can be attributed to a particular design feature of the wuig 
under consideration. Generally speaking, at low Reynolds numbers and no sweep, the laminar boundary layer will 
separate under adverse pressure gradient, the subsequent expansion of the flow will reattach the boundary layer as 
turbulent, and a laminar separation bubble (Fig. 4) will be contained in-between the two regimes. ’ This bubble can 
be identified on thermograms as a narrow, elongated region of low heat transfer preceding the turbulent regime. At zero 
to moderate sweep and higher Reynolds number (but outside the critical conditions for attachment line transition), the 
transition will be caused by Tollmien-Schlichting waves characterized on the thermograms by a relatively smooth 
transition line (Fig. 5a). 68 As the sweep angle increases beyond a certain limit, the cross flow mechanism will prevail, 
with the transition front on thermograms resembling a "saw-teeth" pattern (Fig. 5b). 

The main obstacle to getting wider acceptance for this method, especially in the aeronautical industry, is that models are 
usually made of aluminum or stainless steel of high quality surface finish. The problems mentioned with this type of 
models in the previous section, i.e., high thermal diffusivity, and low emittance and high reflectance, only render the 
method more difficult to apply. In some cases, the problem may be solved satisfactorily by applying an insulating film, 
1 or 2 mm in thickness, in specially machined grooves. 69 Crowder 66 suggested as a partial solution to this problem to 
cover the model with an insulating paint, and to spray liquid nitrogen into the tunnel, upstream the test section. The 
rapidly cooling flow produces on the metallic surface a thermal effect that is strong enough to make transition visible 
(Fig. 6). 

Viewed in the perspective of a very few years, it seems that the IR imaging technique opened new horizons in the 
laminar flow research and applications. The method was used extensively in wind tunnel and flight test programs in 
Germany by the DLR and Deutsche Airbus, 68,70 " 74 in France at ONERA, 6 ,75 and in USA at NASA, 55 ’ 67 to investigate 
the occurrence of the various transition modes. The DLR / Deutsche Airbus program was a particular large-scale effort. 
It combined flight tests on small- and medium-sized airplanes with laminar flow gloves mounted on the wings, with 
wind-tunnel tests performed in different facilities. As part of this effort, an extensively instrumented VFW-614 twin jet 
medium-sized aircraft designated the Advanced Technologies Testing Aircraft System (ATTAS) was dedicated to laminar 
flow research over a period of a few years (Fig. 7a&b). The unique capability of IR imaging to provide repetitively 
global views of the transition pattern on the wing, allowed an immediate identification of the prevailing mechanism 
under given test conditions (Fig 7c). The combined use of IR imaging with hot films throughout the program helped 
elucidate the correlation between the actual mechanisms of transition development and numerical predictions designed to 
verify the e N method. At NASA, the effort focused on comparative flight testing of different competing, or 
suplementing, methods for transition research. 55 At ONERA, the technique was used to confirm designs specifically 
aimed at obtaining significant runs of laminar flow on aerodynamic surfaces. 69 Unfortunately, none of these experiments 
illuminated the question surrounding the identification of the precise location in the transition process that is detected by 
this technique. 

The main drawback with transition studies made in conventional wind tunnels is the scale of the models and the resulting 
inconsistency between the Reynolds number of the model and that of the real flying airplane. This discrepancy is pivotal 
in the search for ground testing methods capable of generating high Reynolds number flows. There are two main reasons 
for concern about this subject. The first is that transition, being a viscous effect, is critically influenced by the Reynolds 


8 



number of the airplane in flight. The second is that as of yet, there are no proven scaling methods for flows especially 

iTTter" 6 For these wind ’tunnels, non-intrus.ve, global and productive transition detection tec^ques are of 
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7. Heat Transfer Studies 


Relative to the obvious qualifications of the method for heat transfer studies, the scarcity of applications is relatively 
surprising. 

Meroney made a wind tunnel study of convective heat transfer of buildings using an actively heated model. The heat 
transfer coefficients were estimated, but lack of a data acquisition and processing system hampered progress on the 
work. Page et al. 88 investigated the flow of a radial jet on a stagnation surface. Using a heated jet impinging on a cold 
surface, they inferred the heat transfer occurring in the stagnation region from thermograms of the heated surface. 
Carlomagno and de Luca, 19 * 89 performed a similar investigation with interest in the heat transfer of both a single jet, 
and an array of jets. In this case, however, a cold air flow impinged on a heated metallic foil. Another jet impingement 
experiment was performed by Eppich and Kreatsoulas. 90 In their case, the surface upon which the jet impinged was a 
thin layer of an insulating material laid atop a constant temperature substrate. Therefore, they could deduce the heat flux 
across the top layer directly from the measurements of the surface temperature, the thickness of the top layer and the 
constant temperature of the substrate. Gartenberg and Roberts 9 proposed a method of mapping velocities by measuring 
the temperature distribution along a thin electrically heated wire placed across the flow of interest. The velocities can 
thus be derived from established, convective heat transfer correlations, provided the error measurement of the IR imaging 
system is known, especially when the wire displays high temperature- gradients in the scanner field of view, Spence 9 ^ 
examined the possibility of determining heat transfer coefficients on vanes used for missile thrust vector control from 
wind-tunnel tests on a model vane under identical flow conditions. Henry and Guffond 92 made temperature 
measurements on a helicopter blade in an icing tunnel to validate computer code predictions of de-icing effectiveness of 
blade heaters. Other examples of heat transfer applications include convective heat transfer from a heated cylinder at 
moderate Reynolds numbers, 93 and external mapping of the temperatures of a gun barrel during firing. 94 


8. A View to the Future 


A review of past achievements must end with a view to the future. Further progress in this field will be contingent on 
two factors: the technical developments in IR imaging technology, and the continuation of experimental exploration and 
aerodynamic research based on this technique. 

Current trends in the electro-optic and semiconductor industries point to future IR imaging systems incorporating focal 
plane detector arrays. Such systems will offer improved performance in terms of spatial resolution, temperature 
sensitivity, and frame rate generation. Optionally, they will trade higher frame rates for lower temperature sensitivity. 
In parallel, as increasingly powerful personal computers become available, they will be used to increase the frame rate 
acquisition and storage. This option will allow capture of faster thermal transients and promote the application of the 
method in very short test duration facilities. 

The testing needs at high Reynolds numbers in cryogenic wind tunnels at low temperatures will require imaging systems 
built around detectors sensitive to longer wavelength IR radiation, up to 30 micrometers. 95 In parallel, there will be 
requirements for window materials with high transmittance in selected bandwidths; good mechanical, thermal, and 
chemical properties; and compatible antireflective coatings. 

At the hot-end of aerothermodynamic research, nonequilibrium, reactive, and high-temperature flow research will require 
multi-spectral IR imaging for measurement of wall thermal response, flow visualization, and tracing of the evolving 
concentration of chemical species through their absorption bands. 
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In «odv».m.c tenting, the ehnllenge is either to solve the tmnsition delation on meullic models or to comp om.se w, h 
othm' materials, either «. inserts .. locrtions of interest, o, as model skins. The exact kcahon where the hansttton , is 
indicated on themtogtams is still far from being defined. More comparative work is needed to define which phase of the 
transition process is detected, and to determine what flow, substrate, and !R imaging system parameters are that in 
the identification of that location. 

In general each flow feature has a thermal signature of its own that should be identifiable on thermograms of models of 
interest As the use of this technique expands, new features appearing on thermograms will puzzle engineer in their 
auest to identify their cause. Separation, shock waves, and vortices have already been observed and identified in the 
past, and continuing research will lead the 1R imaging technique to become a widely accepted diagnostic tool in 

aerodynamic research. 


9. Conclusions 

In the years that have passed since 1967, the infrared imaging of aerodynamic surfaces gained recognition as an 
experimental tool with umque capability for fast mapping of surface temperatures. This capability is used in hypersonic 
research to determine hea^ flux distributions from local temperature measurements and in subsonic research to give 
i r^vit of 7e boundary layer transition to turbulence. With time, other flow features were identified on 
such as shock waves, separation, reattachment, and vortices. The use of the technique was gradually 
expanded to other disciplines, most notably propulsion and basic heat transfer research. 

mi , , - mie had a relatively slow start, due to sporadic experiments, obscure documentation, and lack of data 
processing capability. When viewed in the perspective of the long period of time it took to mature and get recognition, it 
Lems the^ technique was ’invented" too soon. However, the ever-increasing number of publications reporting useoflR 
imaging in aerodynamic research is a convincing proof that the method has passed through the development stage and 
assuming the status of a routine experimental procedure. 
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Fig. 1. a) Space Shuttle forebody model in Tunnel 
B at AEDC von Karman gas dynamics facility 
during hypersonic heating test, b) Typical heating 
pattern on the windward side as captured by the 
IR imaging system. 5 

Courtesy NASA, reprinted with permission. 



Fig. 2. Space Shuttle Columbia atop the Mobile 
Launcher Platform. The Shuttle Infrared Leeside 
Temperature Sensing (SILTS) pod is installed 
atop the vertical stabilizer. 

Courtesy NASA, reprinted with permission. 
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Fig. 3. SILTS data of Space Shuttle heating 

o p rr or t;;i g PI :r^' ,c re ' ent^v • projec,ed r:VM Chara f" s 'ic ,ra “ ,no " 

Courtesy D. Throckmorton and V. Zoby, NASA flow^CFH^ nart a* 1 " 8 waves , b) cross 

Langley, reprinted with permission. ’ anTliSf area T* f in ^ ic . ates laminar re 8 ime 

and light area indicates turbulent regime. Notice 

the saw-teeth pattern characteristic to CF 
instability. 

Courtesy G. Gauffre and V. Schmitt, ONER A 
Chatillon, reprinted with permission. 



Fig. 4. Transition via laminar separation bubble. 

The bubble is indicated by the colder (darker) Fig. 6. Turbulent wedge and transition captured 
straight strip between the laminar (dark) and on a painted stainless steel model by injecting 
turbulent (light) areas. liquid nitrogen into the wind-tunnel flow. 6 ° 

Courtesy A. Quast, DLR Braunschweig, reprinted Courtesy J. Crowder, Boeing, reprinted with 
with permission. permission. 
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Fig. 7. The VFW-614 Advanced Technologies 
Testing Aircraft System, a) Schematic drawing; 
the laminar flow glove, marked gray, is outboard 
the engine pylon.' b) Close-up view of the 
wing area; notice the laminar glove and the IR 
cameras locations. 1 c) In-flight transition 
pattern induced by Tollmien-Schlichting waves 
and insect contamination (turbulent wedges). 
Courtesy A. Quast, DLR Braunschweig, reprinted 
with permission. 
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temperature profile vs. square— root of distance from the 
leading edge. Re x = 1.233k 10^ X(m) 

Data points of three different experiments. 
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The model, assuming a general heat flux distribution chordwi 
ical solution, was formulated by Klein and Iribus (1952). 
to the constant heat flux assumption, the solution reads 


- 1/3 f 'a 1/2 1/2 


0 (x) = 2.2019 ^ — Pr ij]— 

a 


= const • x 


1/2 


/ G 1' 


se ana 

v^i h u ' i 


where 9(x) is the difference between the local wall temperature and tn c 
freestream air static temperature, q“ is the wall heat flux, k a is tne ai r 
conductivity, Pr is the air Prandtl number, ^ is the air kinematic viscosity, 

U is the freestream air velocity and x is the distance from the leading edge. 
According to this model, when the flat plate temperature distribution 
corrected by the ambient air temperature is plotted against the square-root of 
the distance from the leading edge, the slope of the line is an indication of 
the heat flux convected into the airstream. When this plotting is done for 
the present case, as shown in Fig. 5.5, its analysis shows that although the 
data is well correlated ( r 2 is 0.945 with a variance of 0.38) the assumed 
linear behavior of the data is not uniform all along the coordinate scale. A 
closer look suggests one linear regression for the first 14 stations and 
another for the last 18 stations, leaving 6 stations in between as a 
transition zone. Following this approach, the variance of the data goes down, 
for the first group to 0.1 and for the second group to 0.05, which is 
significantly better. The question is what causes this change in the pattern 
behavior. An answer lies in the actual physical process taking place with 
this experimental set-up and its departure from the assumed model (see Eq. 
5.1). As shown in Fig. 5.6, the theory assumes a flat plate "pumping" a 

constant amount of heat flux into the airstream. For a desired constant value 
of this heat flux, a certain prescribed temperature distribution should be 
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a temperature rise proportional to the square root of the distance from 
the leading edge. 

(3) The leading edge experiences an abrupt increase in its temperature, due 
to its sharpness. This geometrical feature causes a finite amount of 
heat to be absorbed by a theoretically infinitesimal substrate mass, with 
the consequent real increase in temperature. The effect is further 
augumented by the fact that the heating wires being wrapped around the 
plate, the leading edge is well heated from all around its contour. 

(4) The temperature distribution exhibits some irregularity around the 
coordinate x=8 cm. This is the region where the thermocouple junction 
was placed, a fact that caused a slight "bumpiness" (about 0.1 mm high 
and 8 mm long streamwise), on the overlaying duct tape surface. This 
geometrical feature initiated a directional emittance factor, which is 
usually lower than its normal counterpart, especially in the range of 
local inclination between 60° to 90° (Sparrow and Cess, 1978, p. 55). 
This effect, not being taken into account by the system's software, is 
ultimately (and falsely) interpreted as a local drop in temperature. 

As a result of these considerations, it was decided to proceed with the 
analysis of the experiments, using only the data points starting at the 
coordinate x = 1.9 cm from the leading edge (the first 1.52 cm from the 
leading edge were ignored). The numerical values of these coordinates 
resulted from the fact that the data was extracted at each third consecutive 
pixel, which on the plate is equivalent to 3.8 mm. 

In search for a model against which the experimental results can be 
evaluated, the physical reality suggests using the flat plate laminar boundary 
layer model with constant heat flux at its surface as the prescribed boundary 
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Fig. 5.4 Constant power heated Hot plate experiment 
DISCO II® output with longitudinal and 
latitudinal temperature distributions. Air 
flows from left to right 
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results ensured that no further benefit would have come out of analyzing more 
profiles, except perhaps to establish a statistical variation measure for the 

data set. 

5.3 Results and Discussion 

The chordwise temperature profiles as obtained from the last three runs 
are shown in Fig. 5.3. The data in Fig. 5.3 was extracted from computer 

stored thermographies, an example of which is shown in Fig. 5.4. Note the 

longitudinal (chordwise) temperature distribution under the thermograpny and 
the latitudinal (spanwise) temperature distribution at three dif erent 
stations along the chord, to the right of the thermography. The wind-tunnel 
air temperature for runs 1 and 2 was 28.6°C and for run 3, was 29.0 C. The 

actual wind tunnel air velocity was 20.9 m/sec during all three runs. The 
camera's field of view in this configuration was 16.5x16.5 cm, and the 

Reynolds number of that area of the plate scanned by the camera varied between 
zero (at the leading edge) to about 217000. This means that the phenomena 
addressed should be all under the laminar boundary layer regime, Re crit = 
500,000. 

The data as presented in Fig. 5.3 displays the temperature profile along 
the plate reduced by the airflow temperature, versus the chordwise 
coordinate. A first evaluation of these profiles reveals the following: 

(1) The information gathered from all three runs was quite reproducible, and 
the data points obtained from the different experiments almost coincide. 

(2) The general qualitative behavior of the surface temperature distribution 
meets the expectations based on the laminar boundary layer theory, i.e., 
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on the heat transfer process from the flat plate to the fluid and in general 
is disruptive to infrared measurements. The other concern was that long 
periods of heating, would affect the substrate to an unknown degree resulting 
in long cooling times between runs and i rreproducible results. In this 
respect, the thermocouple used for emittance calibrations proved to be of 

significant help, in establishing the substrate initial conditions for each 
test. 

According to the above mentioned guidelines, a typical experiment would 
run for less than one minute with the time split between the following steps. 

(1) Wind-tunnel starting and achieving steady-state operation (approx. 15 
secs. ) 

(2) Turmng-on the power supply, and recording of the data. Turning-off the 
power supply (approx. 20 secs.) 

(3) Keeping the wind-tunnel running to cool the plate (approx. 15 secs.) and 
stop its operation. 

The working point of the various systems participating in this experiment 

was determined in previous runs. The nominal test conditions were established 
to be: 

(1) Wind tunnel air speed: 21 m/sec 

(2) Electric current through the heating wire: 345 mA. Total resistance of 
the heating wire: 135 Q 

A total of 20 tests were performed. Of these, eleven tests were 
recorded. Since the results of all the runs were very much alike, only the 
last three were picked-up for data analysis. The reproducibility of the 
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5.2 Cutaway ol the heated Hat plate 



stability principle, 
within 0.5°. 


The accuracy of both devices 


was established to be 
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The target flat plate was placed in the wind-tunnel at +0.5 degree angle 
of attack in order to be on the safe side regarding transition. The active 
heating of the target plate was achieved locally, at one third of the span, by 
wrapping chordwise a 0.00397" Constantan wire three times around the plate at 
2 mm pitch, and connecting it to a power supply. To ensure uniform surface 
texture, roughness and heating, these wires and the entire surrounding area in 
the field of view of the camera were covered chordwise with three adjacent 
strips of 2" wide duct tape. At 85 mm. from the leading edge, a Chromel wire 
was soldered to the Constantan wire and its lead was taken straight to the 
back of the plate through a very small hole that was drilled in the wood. As 
shown in Fig. 5.2, this layout permitted the experimenter to carry-out 
emittance calibrations of the duct-tape surface by providing the true surface 
temperature. Hence, the Constantan wire played a dual role in this experi- 
ment, either as a heating element when connected to the power supply, or as a 
thermocouple element when connected to a digital thermometer. The desired 
function was selected by switching between the two circuits. 

5.2 The Experiment Description 

The guiding rule in performing these experiments was to make them as 
short as possible. There were two reasons for this; one, concerning the wind- 
tunnel operation, the other concerning the substrate participation in the heat 
transfer process. Since the wind-tunnel used in this experiment is of closed 
circuit type, the longer the operation period the higher is its air 
temperature. At room temperature, the initial rise of the air temperature in 
the wind-tunnel is about 1.2°C/min. This phenomenon has 


& negative influence 
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Fig. 5.1 Layout for IR surface measurements of a 
heated flat plate in a uniform freestream 
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Chapter 5 

THE FLAT PLATE EXPERIMENT 

Following the heated wire experiments that essentially involved only one 
mode of heat transfer namely convection, the attention turned to a more 
general type of experiment in which all three modes were involved. The choice 
was a laminar flat plate experiment which was designed to measure surface 
temperatures on an actively heated surface subjected to a uniform free-strean. 
flow. The purpose was to produce an "aerodynamic" target whose temperature 
distribution as measured by the infrared imaging system could be assessed by 
analytic results from velocity and thermal boundary layer analysis. 

5.1 The Experimental Set-Up 

The target flat plate measured 32 cm. chordwise, 90 cm. spanwise, was 
made of wood, and had a sharp leading edge of the type assumed in the Blasius 
analysis of the boundary layer. The IR imaging camera was placed laterally at 
55 cm. from the target, behind a second plate (parallel to the first one!, 
with a circular hole through which the camera could view the target. This was 
done in order to prevent the vortices shedded from the camera from interfering 
with the flow over the target plate. The experiments were carried out in a 
0.9x1. 2 m ( 3 ' x4 ' ) low speed wind-tunnel, where the entire assembly was placed 
vertically as shown in Fig. 5.1. The zero angle of attack of the plates was 
determined by using a cylindrical two-hole differential pressure probe 
connoted to a U-shape glass manometer. In order to eliminate errors 

associated with geometrical misalignments, the zero angle of attack was 
established from the mean of two readings, with the pressure probe at its 
normal position and then rotated 180° around its axis. Separately this 
measurement was done with a vane-type indicator working on the weathercock 
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in the wider context of better understanding the complexity of convective heat 
transfer processes. 
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velocity behavior u to be quite sensitive to the scatter of the temperature 
measurements. In spite of this fact the results of the experiments with the 
higher heating rates (and thus the best correlated data) prove that the 
concept of deducing velocity information is viable. 

The encouraging aspect of this approach is that as the Reynolds number 

increases with the velocity, so does the heat transfer coefficient, on a power 

law basis (Morgan, 1975). The net result of this feature is -o increase 

substantially the accuracy of the velocity predictions from the temperature 
measurements. 
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This set of experiments showed that using an infrared imaging system, one 
can detect local temperature variations, and hence Nusselt number variations, 
that may occur over a very limited portion of the camera's field of view. The 
concept proved to be applicable in the core region of the flow as well as in 
the entrained part of it. The capability of accurately tracking temperature 
variations was found to be in this case as high as 350 K/m. A question still 
to be answered is whether or not this limit is set by the actual value of the 
temperature gradient, or by the variation of the temperature in the field of 
vie. of the IR imaging camera. The inverse approach of using temperature 

measurements in conjunction with heat transfer correlations to deduce local 
velocities was shown to be possible in principle. This method shows promise 
for use in aerodynamic research: (1) provided the wire overheat is high 

enough to get good signal to noise ratio measurements, (2) the temperature 
gradient problem is solved by appropriate calibrations, and, (3) more accurate 
forced convective heat transfer correlation may be found. The last condition 
should be viewed not only in the perspective of this particular research, but 
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decreases, the discrepancy between predictions based on Morgan's correlations 
and the data obtained increases. It may reach as much as 30% for Ra numbers 
as low as 300. In the present case, the Ra numbers were even lower, typically 


The deduction of the velocity profile from the measured temperature 
profile through Nu-Re correlations is also possible. This time, the 
process is reversed and assuming the correct correlation is available, the 
known differential temperature is used to calculate the velocity that would 
produce it. This approach is schematically shown in Fig. 4.9 and quantita- 
tively described in Appendix A2 item 5. The results for the first two 
experiments (lower heating rates) are presented in Fig. 4.10 while the results 
for the last two experiments (with the highest heating rate) are presented in 
Figs. 4.11 and 4.12. For the predicted velocity values, Eq. (4.3) was 
reformulated to express the local jet velocity U(r) in terms of the local 
temperature difference A T(r) between the heated wire and the air 


0.909 v 


U(r) = 


[■ 


it k AT(r) 
d 


2.604 

] 


(4.4) 


Once again the reader is reminded that the two most peripheral points in each 
graph should be regarded more as a mathematical feature (arising from the 
assumed parabolic distribution) than an actual physical reality. A comparison 
between Figs. 4.10-4.12 only emphasizes how important is in this case to get 
well correlated temperature measurements before proceedings with the data 
processing (see Table 4.2). The explanation of the higher scatter of the 
deduced velocities (Fig. 4.10) than that of the temperature data from which 
they were deduced (Fig. 4.6) lies in the behavior of the Nu-Re correlations. 
These correlations assume a behavior like, Nu ~ Re n 0 < n < 1. In the former 
case one gets AT~U“ n while in the latter U~(AT)"^ n , thus causing the deduced 
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Table 4.4 Mean temperature difference between the heated wire 
under free convection and the ambient air as a 
function of the wire heat dissipation. 


II 

, % 

(w/m wire) 

AT predicted (K) 
Nu-Gr corr.* 

AT (K) 

Experimenta 1 

3.47 

65 

25.6 

4.71 

36 

34.5 

6.17 

110 

44.8 

9.33 

161 

63.9 


♦See Appendix A2 item 3. 
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processes on the raw (detector) signal, that practically nothing is known 
about the transfer function operating on it. First, there is the IR imaging 
camera, then its accompanying analog processor and finally the computer 
digitizer and its software. (Several attempts to get some information about 
these transfer functions from the suppliers were not answered). Therefore, 
there is no choice left but to look at this system as a "black box" and to 
assess its performance by controlled experiments, even though this task is 
double complicated because of the uncertainty induced at the same time by both 
the measurement apparatus and the reference data. In any case, the trend in 
Fig. 4.8 shows that the accuracy of the deduced Nusselt numbers from tempera- 
ture measurements performed with the infrared imaging system is adversely 
influenced by the existence of a temperature gradient at the location of 
interest. It also confirms that this system may be able to track accurately 
temperature variations up to about 350 K/m which is in good agreement with the 
findings of Stallings and Whetsel (1984). Furthermore, according to the 
evidence presented, at high temperature gradients the indicated temperature is 
always lower than the actual one. This suggests that it is possible to 
produce a calibration curve similar to Fig. 4.8 and use it to correct the 
experimental measurements prior to processing the data. 

The uncertainty with respect to the accuracy of heat transfer correla- 
tions is further emphasized by the results measured on the uncooled portions 
of the wire. Although the free convection heat transfer was out of the scope 
of the present study, it was interesting to observe that on those segments, 
the predicted Nusselt number as regressed by Morgan, grossly underestimates 
the actual heat transfer process, as it is shown in Table 4.4. 

The verity of the present findings are indirectly supported by Fand et 
al. (1977) where it is shown that as the Rayleigh number (Ra ■ Gr Pr) 
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heating rate is known and constant, the dissipation is effected solely by 
convection, and the actual wire temperatures are higher than the measured 
ones. Therefore, Nusselt numbers that are deduced based on this kind of 
information should be corrected to actual values by substrating the error 

inherently associated with this type of measurements. 

One question that remains open on this issue is whether the infrared 
camera limitation to track high temperature gradients is a function of the 
absolute value of the gradient itself (K/m), or of the variation of the 

temperature in the instantaneous field of view (K/IFOV or K/rad). To answer 

this problem, it would be necessary to undertake another set of experiments, 
deviating from the objectives set-out for this research. However, the 
resemblance between Figs. 4.8 and 3.2 relating the minimum detectable 

temperature difference is remarkable. Indeed, the prediction made on the 

basis of the IR imaging camera operational characteristics about its inability 
to accurately resolve high temperature gradients turned out to be true. 

The quantitative estimation of the performance of the infrared imaging 
system in this case is a very difficult task. The first difficulty refers to 
the reliability of the heat transfer data chosen to be the reference. As 

Morgan points out, "The percent coefficient of variation of the Nusselt number 
for a given Reynolds number for the experimental data varies from 10% to 29% 
depending on the Reynolds number, whereas that for the various correlations 
varies from 10% to 46%". Based on the present measurements, and referring 
only to those performed in near zero temperature gradient regions, it was 
found out that the experimental Nusselt number based on infrared imaging 
system measurements could further vary by as much as 3.3% from the predicted 
one by Eq. (4.1) to which Morgan’s remark refers. The second difficulty 
refers to the fact that the experimental data set resulted from a series of 
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In this case, it will be easier to evaluate the influence of the 

temperature field gradient on the measurement accuracy, because of the 

position dependence of the temperature gradient along the wire. For this, it 
was assumed that the actual wire temperature distributions could be 

approximated by a cosine-type function between X = 0 at the maximum wire 
temperature and X = x/2 at the minimum wire temperature (0 and .radians 
respectively). This approximation function was suggested by the experimental 
data shape. While doing this, it was assumed that the temperature 

measurements are essentially correct in regions with no sharp gradients 
(Stallings and Whetsel, 1982). The approximating functions are presented in 

Appendix A3. The final result is presented in Fig. 4.8, where the estimated 

error in the measured Nusselt number relative to the expected values from Eq. 
(4.1) are shown versus the temperature gradient along the wire at the location 
of measurement. The Nusselt number error shown in Fig. 4.8 was uniformly 

adjusted by the difference between the experimental and the correlation values 
found at the point with zero temperature gradient. This small discrepancy ' 
(see Fig. 4.7) was assumed to be a systematic error applying all along the set 
of data belonging to the same experiment. The data points for Fig. 4.8 were 
taken from the last two experiments (with the higher heating rates) were the 
signal to noise ratio was previously found to be high enough to produce well 
correlated data. Oata points for the two most peripheral points in each 
experiment, to which reference was made previously, have been excluded. The 
line in Fig. 4.8 should not be interpreted as the regression curve on the data 
points. It is rather shown there with the purpose of indicating the trend of 
the findings. The conclusion of this analysis is that Nusselt numbers values 


deduced from the IR imaging system temperature measurements will always be 
higher than the actual ones. The explanation lays in the fact that the 
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Table 4.3 Nusselt number comparison for the laminar flow jet 
using the heated wire technique for experiment *4 
(highest overheat). 


Pi xel 
# 

Re* 

[Actual ] 

Nu( correlation) 
[= 0.795 Re 0 * 384 ] 

Nu 

[Actual] 

Nu(act)-Nu (cor) 

Nu(cor) ‘' 0j 

0 

17.013 

2.361 

2.384 

0.9 

+1 

16.508 

2.333 

2.357 

1.0 

-1 

16.508 

2.333 

2.337 

0.2 

+2 

14.976 

2.248 

2.312 

2.9 

-2 

14.976 

2.248 

2.337 

4.0 

+3 

12.423 

2.092 

2.268 

8.0 

-3 

12.423 

2.092 

2.244 

7.3 

+4 

8.892 

1.840 

2.180 

18.4 

-4 

8.892 

1.840 

2.163 

17.6 

+5 

4.297 

1.397 

2.029 

45.2 

-5 

4.297 

1.397 

2.009 

43.8 


^Decreasing Reynolds number coincides with increasing temperature 
gradient due to the velocity distribution of the flowfield. (See 
Figure 4.7 and Appendix A3) 
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influence on the results. The radiation effect was found to be negligible. 
Because of the moderate wire overheat values, the emittance was assumed 
constant and, according to Kirchhoff rule, the absorbtance could be taken to 
be the same as the emittance (Eckert and Drake, 1972, p. 612). Thus, even for 
an overheat of, say, 80K above the room ambient the radiation effect turned 
out to be three orders of magnitudes lower than the generated heat. In order 
to estimate the possible heat conduction contribution, the measured tempera- 
ture profile of the experiment with the highest overheat was approximated by a 
polynomial function that indicated a maximum temperature gradient of 2524K/m. 
With a mean wire thermal conductivity of 50 w/mK, the conductive loss for a 
wire element between two adjacent sampling points is still, two order of 

magnitudes lower than the generated heat and the influence of this mechanism 
can be discarded too. 

Table 4.3 presents in a tabular form the departure of the Nusselt number 
as deduced from the temperature and heating rate data of experiment #4 from 
the Nusselt number as predicted by Eq. (4.1), and assuming the parabolic 
velocity distribution inside the jet. The positive departure values exhibited 
by the results raises the suspicion that the infrared imaging system may have 
difficulty in accurately mapping high gradient temperature fields. Stallings 
and Whetsel (1982) came to a similar conclusion using a similar camera 
(Thermovi s i ort* 680) to map the temperature variation along a bar with an 
almost constant temperature gradient. Their data show that while the system 
could track temperature gradients of 109K/m (5R/in) with a ±1% error, at 
262K/m ( 12R/ in) this error increased to -1.5%, but they did not pursue this 
issue to determine the dependency of the error on the temperature gradient of 
the measured field, as may be done in the present circumstances. 
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Table 4.2 The coefficient of determination r^ for the temperature 
measurements made with the infrared imaging system using 
the heated wire technique at various levels of heating 
rates, exposed to a circular jet. 


Experiment 

0 

Heating rate 
[w/m wire] 

★ 

r 2 

rot i 

1 

3.47 

41.0 

2 

4.71 

53.1 

3 

6.17 

80.1 

4 

9.33 

98.2 


♦Adjusted for degree of freedom. 
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where U(r) is the local air velocity inside the jet, U max is the maximum air 
velocity at the center of the pipe, r is the radial coordinate and R is the 
internal radius of the pipe. In this case, 'J max = 4 m/sec. The final result 
that correlates the temperature difference between the heated wire £.T(r) and 
the airflow (see Appendix A, Eq. A7) is given by 


AT(r) 


0.7 95 7i k. 


20 


max 


0.384 



(4.3) 


where q is the electrical heat generation per unit length of wire, k a is the 
heat conductivity of the air and v a is the kinematic viscosity coefficient 
of the air. The Reynolds number at the computed points varied roughly between 
4 to 20. Figure 4.6 indicates that increasing the wire overheat value 
improves the signal to noise ratio and the resulting acquired data is better 
correlated. The quantitative evidence is presented in Table 4.2 where r is 
the (statistical) coefficient of determination (see Appendix A, Eq. A. 4). On 
the negative side, high overheating may cause changes in surface emittance, 
excessive wire elongation (with resulting vibrations) and thermal contamina- 
tion of the flow. Thus, the overheat value should not be set to exaggerated 
values. A second observation relates to the fact that the measured temperature 
values systematically depart from the expected values according to Eq. (4.3.). 
These trends become more evident in the data from the experiment with the 
highest heating level presented in Fig. 4.7, where the wire temperature 
difference is graphed in order to make the results more meaningful for heat 
transfer evaluations. 

Since the correlation given by Eq. (4.1) accounts only for the forced 
convective heat transfer, the magnitudes of the other two heat transfer modes, 
namely conduction and radiation, were estimated in order to assess their 
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cooled section of the wire length indicating the jet width will be the same 
too, under the assumption of negligible radiation and conduction. Hence, the 
"wavelength" of the wire exhibiting a varying temperature distribution will be 
the same in all four experiments. The validity of this assumption is examined 
next . 


4.4 Results and Discussion 

The raw data of the initial experiments at the three lower heating rates 
are presented in Fig. 4.6. The darkened symbols are given for comparison and 
represent the expected temperature distribution along the section of the wire 
directly exposed to the laminar flow. The respective temperatures were 
predicted from the relation, 

Nu = 0.795 Re 0,384 . 1 < Re < 35 (4.1) 

which is a cross-flow correlation deduced by Morgan (1975) through statistical 
analysis of some 200 independent experimental and theoretical studies. 
Although Eq. (4.1) refers to two-dimensional flows, and the circular jet is a 
three-dimensional flow, it is still possible to use this equation in the 
present application because the wire diameter is almost three orders of 
magnitude thinner than the jet. Thus, it is assumed that the wire is on the 
plane of symmetry of the jet, where the flow, locally can be assumed two- 
dimensional. The air velocities to be used in Eq. (4.1) were calculated 
assuming parabolic velocity distribution as predicted by the viscous-flow 
theory 

u l r > - u ,nax t 1 - (£> 2 1 


(4.2) 



Fig. 4.5 Comparison of the normalized velocity distribution 
in a circular turbulent jet with a fully developed 
laminar flow in a pipe (parabolic distribution). 
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the periphery, the jet will experience an enlargement due to the mixing 
effect. The magnitude of this singular enlargement cannot be predicted in a 
straightforward way from the theory. Due to mathematical boundary conditions 
required for solution, it may be assumed that the jet issues from a round slit 
in an infinite wall. One of the techniques to match the theory with the 
reality is to assume the nozzle moved back with the wall a virtual distance 
that is calculated through a suitable matching criterion. For example, it can 
be done by matching the volumetric rate of flow between the actual jet and the 
assumed one (Blevin, 1984, p. 232). This analysis indicates that the actual 
flowfield behavior around the wire lies somewhere between the parabolic and 
the circular jet velocity profiles. The two normalized velocity profiles, 
parabolic versus circular jet, are shown for comparison in Fig. 4.5, where the 
positions of the indicated points on the theoretical parabolic distribution 
coincide with the actual temperature measurement locations as shown in Figs. 
4.6 and 4.7 with an error smaller than 2.7% (These figures will be discussed 
in detail in the next section). This comparison suggests that the uncertainty 
with regard to the actual flow velocity may be too large at the two 
(synsnetrical ) most peripheral points and th^ data there should be treated 
cautiously. Another point emerging from Fig. 4.5 is that the effective radius 
of the circular jet flow is about 1.8 times the radius of the pipe. This fact 
is directly supported by the experimental evidence presented in the next 
section in Figs. 4.6 and 4.7, where the cooled section of the wire is larger 
than the internal radius of the pipe. Assuming that the extent of the cooled 
section is an indication of the extent of the actual flowfield, it can be 
inferred that the capability of the infrared imaging system to map the core as 
well as the entrained flow regions is demonstrated. Indeed, if one assumes 
the flow velocities to be approximately the same in all four experiments, the 
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assumed the same, due to the viscous interaction between the flow and the 
still surrounding air. This interaction causes the well-known entrainment 
effect which affects the flow velocity profile especially at its periphery. 
Since the air velocities near the pipe wall fall to extremely low values tnat 
are difficult to measure, an attempt is made to estimate the velocity profile 
perpendicular to the wire. According to Blevins (1984, p. 231) a circular jet 
is fully laminar for Reynolds numbers up to 1000 and fully turbulent for 
Reynolds numbers above 3000. This Reynolds number is based on the diameter of 
the nozzle from which the jet issues. Regardless of the jet type of flow, 
whether laminar or turbulent, both of them exhibit the same type of behavior 
with regard to the main parameters of interest: the jet wjdtn increases 
linearly with the axial distance (X) while the centerline velocity decreases 
with inverse proportionality to the axial distance (1/X), see Schlichting 
(1968), p. 686. To this, it must be added that in both cases the flow pattern 
fully establishes itself only after a certain distance downstream from the 
pipe exit. In this case of intermediate Reynolds number of 1700, the 
distinction between laminar and turbulent jets is more of a semantic matter. 
The difference is whether one assumes the transition to be an "entrainment" 
phenomenon (laminar) or-a "mixing" one (turbulent). Therefore, and for these 
experiments only, further reference will be made simply to "the circular jet". 
It is assumed that the coreflow is very much like the one in the pipe, the 
fractional change in the centerline velocity at the wire being of the order of 
magnitude of 10 - ^ relative to the nominal velocity in the pipe. This 
assumption was supported by measurements of the coreflow at the mouth of the 
pipe performed with a static Pitot-tube connected to a U-shaped glass 
manometer inclined 85.3° from the vertical. The measured core velocity 
distribution versus the parabolic velocity profile is shown in Fig. 4.4. At 
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correction for its discharge coefficient. The throat diameter of tne nozz'e 
was 0.8 mm (0.0315 ), and the nominal absolute air pressure upstream the 
nozzle was 304510 N/m 2 (31 psig). This pressure was adjusted to account for 
the nozzle discharge coefficient and variations in the air temperature, to get 
a nominal centerline velocity of 4 m/sec, determined from Pitot-tuoe 
measurements (see Sec. 4.3). 


4.2 The Experiment Description 

The experiments were carried-out at four different levels of heating 
corresponding to dissipation r tes of 3.47, 4.71, 6.17 and 9.33 watts/m cf 
Chromel wire. The area of i terest i.e., that part of the wire that was 
directly cooled by the air jet, was about one tenth the field of view of the 
IR imaging camera. This design was made purposely to evaluate the system's 
capability to resolve localized temperature changes in the field of view. 

For each separate experiment, nine consecutive frames were taken at a 
rate of 0.7 Hz and their average was stored as a single frame on tne 
computer’s diskette for later analysis. During each experiment, the range and 
level of the camera were carefully set according to the target luminosity in 
order to minimize the background noise, and at the same time to prevent the 
detector from being under or overexposed (saturated) to the target photons. 
In either case, the output ceases to be proportional to the fourth power of 
the target temperature. 


4.3 The Flowfield 

As previously explained, the wire was placed at the mouth of a pipe 
designed to produce a laminar flow, implying a parabolic velocity profile 
inside the pipe. However, the behavior at the exit of the pipe can hardly be 
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The Chromel wire length was 40 cm and the Constantan wire length was 10 
cm. The thermocouple assembly was hung between two vertical bars 50 cm apart. 
The length to diameter ratio of the Chromel wire was 5250:1, so that the 
influence of the heat conduction (on the central area of interest) to the 
supports can be neglected. The wire was heated by a O.C. electric current 
passing through it. A cardboard screen was placed behind the wire to provide 

a uniform background for the IR imaging camera (see Sec. 2.1). As shown in 

Fig. 4.2, the pipe exit was placed at 2 mm distance from the heated wire, at 

the center of the Chromel section. The air-supply pipe bad a 13.3 mm (u.545 ) 
internal diameter, a 15.9 mm (0.625") outer diameter and a straight section of 
80 cm, thus getting a length to diameter ratio of 60. The nominal mean air 
velocity in the pipe was about 2 m/sec, and the corresponding Reynolds number 
based on the pipe diameter was about 1700, which is well under the critical 
Reynolds number of 2300. The length of the straight section was considered 
enough to provide fully developed laminar flow inside the pipe, just before 
its mouth. According to an approximate analysis by Schiller (Schl ichting , 
1968, p. 231) a length to diameter ratio of 49 could be considered satisfac- 
tory for the actual Reynolds number. According to another analysis by 

Langahaar (Kays and Crawford, 1980, p. 67), this design is acceptable within a 
margin of 5% off the ultimate coefficient of friction for fully developed 
flow. Three fine-mesh screens, (each rotated at 45° with respect to the 
others) were placed at the entry of the straight section of the tube in order 
to break down any large vortices or non-uniformities that might exist in the 
flow. The air was supplied by a standard 100 psi low-pressure system (of the 
type that is commonly in use in laboratories and workshops). The air mass- 
flow rate was adjusted using a pressure regulator and by taking total pressure 
and temperature measurements upstream from a sonic nozzle with the suitable 


Table 4.1 Calibration of the Apparent Wire Emittance 
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Fig. 4.3 Set— up tor wire emitfance calibration. 








39 


Calibrations for surface emittance properties are usually pone ted-ous 
processes. As shown in Fig. 4.3, tne calibrations were performed by aiming a 
working heat-gun at the wire segment including the thermocouple junction, sous 
getting a direct reading of the wire temperature. Previously, the heat gun 
position was carefully adjusted in order to get locally a uniform ; fIrt ) 
temperature distribution of the wire in the region adjacent to tne 

thermocouple junction, within t0.2°C. The hot air temperature distribution in 
that region was double checked for uniformity, independently, „itn another 
thermocouple, and the mean hot air temperature was measured with a mercury 

thermometer. The accuracy of each of these devices was +0.2"c. In parallel, 
the wire was thermographed by the AGA system and the temperature of tne 

Chrome 1 wire adjacent to the thermocouple junction was analysed with the Disco 
II® software on the dedicated microcomputer. The emittance input of the wire 
was adjusted until the temperature result from the computation equaled that of 
the thermocouple readout. Since the Disco II® computer program that was used 
to analyte the thermographies could read only emittance values down to two 

decimal digits, recourse was made to the atmospheric transmittance factor to 
finally adjust the emittance value. This result was checked with another 
software package supplied by AGA and run on an HP-41 CV programmable 
calculator with identical results. Thus, for the 0.0762 mm (0.003") diameter 
Chrome! wire at 0.6 m distance from the 1R imaging camera and moderately 
overheated (say around 350K), this apparent emittance has a value of 0.0244, 
as indicated by the HP 41 CV calculator run program. With the Disco II® 
software, the emittance value was set at 0.03 and the transmittance at 0.82, 
to get an equivalent apparent emittance value of 0.0246. The sequence of 

steps for one such wire calibration is given in Table 4.1. 
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Chapter 4 

THE LAMINAR FLOW JET EXPERIMENT 

Following the promising preliminary evaluation, the next step was to 
assess the capability of the IR imaging method to produce relevant 
quantitative results. For this purpose the research was focused on the 

mapping of a laminar flow air jet using the heated wire concept, which is 
shown in Fig. 4.1. The ordered behavior of this flow offered the convenience 

of comparing well understood experimental and theoretical results. The 

measured temperatures were used in conjunction with heat transfer correlations 
in order to deduce the air velocity distribution within the jet. In 

particular, this flowfield offered the opportunity to test the capability of 
the system to map high gradient temperature distributions. 

4.1 The Experimental Set-up 

In this experiment, the infrared imaging camera scanned the temperature 

distribution along a thin, electrically heated Chromel wire, a small part of 
which was placed diametrically at the exit of a pipe from which a laminar flow 
air jet discharged into the atmosphere. The camera was placed at a distance 
of 0.5 m from the wire. At this distance, it has a spatial resolution of 
0.0034 radians and its field of view was 15x15 cm. As shown in Fig. 4.2, the 
Chromel wire, 0.0762 mm in diameter (0.003"), was part of a Chromel -Constantan 
thermocouple assembly that offered the opportunity to perform in-situ calibra- 
tion for the surface emittance value. The knowledge of this value is critical 
when the objective is to determine surface temperatures and air velocities 
from infrared thermographies. 
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3.6 Sunnary 

Use of infrared imaging systems requires the understanding 0 f their 
principles of operations, as well as the characteristics of the specific 
system which is placed in service. These systems do not measure the 
temperature of the object being scanned. Rather they have a voltage output 
which is proportional to the photons number coming from the target and hitting 
the detector. This input may be attenuated by a multitude of factors, among 
them the surface emittance and the atmospheric attenuation. Therefore, in 
extracting quantitative data from system readings, it is the responsabi 1 i ty of 
the experimenter to account for all the possible influences that may have 
influenced the result. Also, it was shown that the IR imaging systems may 
experience difficulties in mapping high temperature gradient fields. 
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so does its angle and the signal modulation, until the detector response 
reaches a plateau characteristic to resolved targets, as shown in Fig. 3.3c. 

With the above background information, the behavior of the curve for the 
minimum resolvable temperature difference, shown in Fig. 3.2 may be 
understood. As the slits frequency increases, the modulation of each 
individual slit decreases, and one needs a higher temperature difference 
between the two plates in order to get a minimum modulation of the signal , 
enough to be detected by the detector. 

From the practical point of view, this behavior suggests that the .R 
imaging camera may have difficulty in tracking high gradient temperature 
fields. Visually, this may be explained in the following way: Let's 
approximate the temperature distribution of a target by a series of constant 
tempecature strips. As the temperature gradient increases, the width of these 
strips decreases and therefore their individual modulation decreases. Thus, 
as the IR imaging camera scans such an object, the error of the measurement 
will increase. In fact, the graph showing the dependency of the error in the 
measurement of a high gradient temperature field will look very much like that 
in Fig. 3.2, with the abscissa representing the gradient of temperature and 
the ordinate representing the mesurement error. 

Before concluding this chapter, it is well to remember that the infrared 
radiation includes diffraction behavior, which is typical to electromagnetic 
phenomena. This means that the signal modulation will not be full even if the 
target is just resolved, since energy from the surroundings will still reach 
the detector. As a rule of thumb, the target area to be scanned for full 
modulation should be at least three times that of the instantaneous field of 


view ( IFOV ) . 
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3.5 The Resolution Problem 

The essence of a successful design for an infrared imaging system lies in 
the answer to the resolution problem. This problem has two aspects: tne 

minimum resolvable temperature difference and the slit response function. 

'he minimum resolvable temperature difference (MRTD) is the result of a 
standardized test desicied to evaluate the camera performance in tnis 
respect. The test target is a plate with four slits of varible widtn which is 
placed in front of another uniform plate. The temperature of both plates is 
independently adjustable. Sy varying the width of the slits and recording the 
respective temperature difference between the plates, where the pattern 
becomes indistinguishable to the camera, one obtains a graph as shown in Fig. 
3.2. The lowest point of the curve is the temperature difference at zero slit 
spatial frequency and represents the system sensitivity. The test is usually 
run at 30 C, and for the present system, this sensitivity is about 0.2°C. A 
second parameter is the highest possible resolution, which can be represented 
by a vertical line that is asymptote to the rising curve. To understand this 
behavior, one must understand the slit response function (SRF) of the 
camera. As it can be seen in Fig. 3.3a, the test target in this case consists 
of a temperature adjustable plate covered by another plate with a slit of a 
variable width. As the slit width is increased, so is its angle as viewed 
from the detector. This angle can be less, equal or greater than the field of 
view of the detector which is called the instantaneous field of view (IFOV). 
Slits that are smaller than the IFOV are called unresolved. When the camera 
scans an object that is unresolved, the radiation falling on the detector 
comes partially from the target and partially from the surroundings. As a 
result, as shown in Fig. 3.3b, the signal modulation from an unresolved slit 
is smaller than one coming from a resolved slit. As the slit width increases, 
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where e is the spectral emissive power of the black body radiation, h is 
Planck's constant, k is Boltzmann's constant, c 0 is the speed of light in 
vacuum, n is the refraction index of the medium, \ is the respective 
wavelength and T is the absolute temperature. When Planck's formula is 
divided by the energy of a single photon h c q /\ = hv we obtain the spectral 
photon emittance for a blackbody 

2 7tc 

N = 2 (3.2) 

n 2 \ 4 [exp(hc o /n\kT)-l] 

From here, one gets directly the formula commonly used to express the 
calibration curve behavior, 

I = A (3.3) 

1 C exp(B/T-T7 

where I is the "thermal value" corresponding to temperature T and A, B, and C 
are constants (AGA, 1984). The thermal value I is measured in isothermal 

units (IU) which is a practical unit of measurement. The relationship between 
the thermal value and the received photon flux is linear. However, the 

relationship between the thermal value and the object temperature is obviously 
non-1 inear. 

Finally, it should be mentioned that the radiation coming from a certain 
spot to the detector has three components: object radiation, reflected 

radiation and atmospheric participation. The camera has no way to distinguish 
between the three, and its output is an "integrated" one. It is indeed the 
duty of the operator to account for the various contributions causing the 
final output signal. As a corollary, it can be said that the measurement of 
temperatures using infrared imaging techniques is not a method that can be 
applied without expert attention to detail. 
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the IR radiation all along its spectrum, the matching of the detector has also 
to fit the "atmospheric windows". This is the term used to denote the 
wavelength of IR radiation for which the atmosphere is transparent. In the 
present case, the sensor is made of Indium Antimonide (InSb), which is 
classified as a semiconductor element with a spectral response in the 1.7 to 
5.8 u range. However, the spectral response of the camera (as a system) is 
purposely reduced to the 3.5 to 5.6 p by special coatings on the objective 
lens. Its operational temperature detection range is between -20°C to 850°C. 
The field of view (FOV) of the objective lens of the camera is 20°. The area 
projected at any one time on the detector is subtended by an angle of 0.0035 
radians, which is also called the instantaneous field of view (IFOV). With a 
scanning rate of 25 fields per seconds, each field having 100 lines and 130 
columns, the scanning rate is 325,000 IFOV/sec, or 3 p sec per IFOV. With a 
time constant of approximately 0.2 p sec (Levinstein 1977), the InSb detector 
matches the camera needs. Since the camera design probably started from the 
detector characteristics, it is clear that they are the bottleneck for the 
upper limit of the scanning rate. 

The InSb is a linear photon counting photovoltaic detector. These 
detectors are sensitive to the number of photons hitting them (rather than 
their energy), once their energy exceeds the energy threshold of the 
semiconductor material. In general, the effect of cooling the detector is to 
decrease this cut-off wavelength. The calibration curves of the IR imaging 
camera can be derived directly from the Planck's formula (Sparrow and Cess 
1978, p.6) 


2ithc‘ 
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n 2 X. 5 [exp(hc 0 /n\kT)-l] 


e bx< T > 3 


(3.1) 
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As can be seen in Fig. 3.1, the scanning mechanism consists of two 

mutually perpendicular prisms. T*e horizontal one rotates at three ros. Its 
optical output is passed to the vertical one which rotates at 30C rps. Thus, 

we get the "lines" and the "columns" of the scanned scene, which is callea a 

field. Each field has 70 active lines (100 all together) and 130 active 
columns. For visual display reasons, the vertical and horizontal motors are 
synchronized in such a way that four fields produce one interlaced frame. 
With a scanning rate of 25 fields per second, 25/4 completely interlaced 

picture frames per second are thus obtained. The output from the vertical 
prism is focused onto the single element semiconductor detector, wnicn is 
located in the wall of a Dewar chamber, and is maintained at 77K by liquid 
nitrogen. The detector produces an electronic output proportional to the 
incoming infrared radiation. After amplification, the signal is sent via a 

cable to the display unit, where absolute or relative temperature measurements 
can be made. 

Throughout this research, the data was aquired only in the "field" mode, 
to take advantage of the higher scanning rate associated with that mode (25 
images per second). More information about the system and the way it operates 
can be found in the AGA Thermovi siorf 8 782 Operating Manual (1984). 

3.4 Physical Aspects of the Camera Operation 

The operational characteristics of an IR imaging camera are largely 
dependent upon the detector characteristics and performance. Each kind of 
detector has its own narrow range of detection in the IR spectrum according to 
its chemical composition. Therefore, it is important to find out a priori 
whether the expected targets do have a radiation component that matches the 
detector characteristics. Secondly, since the atmosphere does not transmit 
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record single images, or an image that is the average of several images, or to 
record a series of up to 16 images at a time at a frequency of 1.5 Hz or 

lower. A key feature of this program is the gradual assignment of eight tones 
of artificial colors (from black through blue, green, red, etc. to white), for 
corresponding gray shades in the original image, thus producing an artificial 
sensation of the temperatures in the scanned field. The main advantage of 
such a system is the cross-hair cursor that can easily be moved on the screen 
and allows the operator to get precise temperature readings at defined 

locations of interest. Although the Discd® 2.0 was written specifically for 
the BMC IF800 computer, and although the latter has a quite obsolete 

architecture, together they make-up a quite versatile system that offers an 
abundance of processing options, from which only the main ones have been 

reported here. 

3.3 The Infrared Imaging Camera 

The infrared imaging camera converts electromagnetic thermal energy 
radiated from the scene contained in its field of view, into electronic video 
signals. These signals are electronically amplified and transmitted to the 
black and white display unit where the signals are further amplified, and the 
resulting image is displayed on the screen. The camera, or the scanner, 

comprises the following main subsystems: 

(a) Electro-optical scanning mechanism 

(b) Infrared detector attached to a liquid nitrogen Dewar flask 


(c) Control electronics and preamplifier. 
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different requirements from an IR system when it is used for nignt vis, on 
imaging versus imaging for temperature measurements, and the inherent cnanges 
to be made in the former in order to realize the latter. A state of the art 
review of infrared detectors is given by Levinstein (1977). Altnough written 
twelve years ago, the article is still relevant as some of the detectors 
described there are still in use. With regard to the quantitative 
determination of surface temperatures using an infrared camera, Hsieh anc 

Ellingson (1977) propose a method that also includes the surface reflection 
problem. 

Before proceeding with the description of the technical aspects of 
infrared imaging systems, it is worthwhile to draw attention to the 
proceedings of Thermosense V (1982), that contains a few papers concerned with 
the practical aspects of temperature measurements using these systems. 

3.2 The Infrared Imaging System 

The infrared imaging system used in this investigation was an AGA 
Thermovision® 782. The basic system consisted of the scanner and a black and 
white display unit that allows direct or relative temperature measurements. 
The output from the scanner can be directly recorded on a modified vidicon 
recorder (VCR) for later playback. This device was of very little use mainly 
because of two limitations: frames could not be addressed individually; and a 

lack of a time base that could be recorded on the tape. The basic system was 
connected through a data link to a microcomputer BMC IF800 dedicated to the 
analysis of the thermographies. A Discd® 2.0 software package is used for the 
image processing and completely occupies a 5 1/4" floppy disk in a read-only 
storage mode. Up to 36 individual images can be stored on a second floppy 
disk in a read and write access storage mode. The user has the option to 
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Chapter 3 

THE INFRARED IMAGING SYSTEM 

A measurement system can be useful only to the extent that the operator 
understands the capabilities and limitations of his device. Although easy to 
use, the infrared imaging systems are very complex electro-optical devices and 
their output is influenced by a myriad of factors besides the target surface 
temperature. Since, in general, aeronautical engineers are not too familiar 
with these systems, it is the purpose of this chapter to review the main 
points of interest related to their performance. 

3.1 Some Useful References 

Infrared imaging systems work on the principles of radiation heat 
transfer. This type of heat transfer is based solely on the fourth pwer 
temperature difference between the body of interest and its surroundings, and 
takes place through electromagnetic radiation at wavelengths between 0.3^ to 
50iu (where n = 1CT 6 m), known as the infrared spectrum. As the name 
implies, the IR imaging systems detect this radiation and through electronic 
processing of the signals, produce an artificial video picture of the area 
scanned in which darker shades are associated with lower temperatures and 

lighter shades with higher temperatures. 

The infrared imaging systems used for surface temperature measurements 

were developed during the sixties as derivatives of the military passive night 
vision systems. Although the requirements from the two types are somewhat 
different, the operational principles are the same. Lloyd (1975) gives an 
excellent overall account of the subject as viewed from his perspective as an 
engineer who worked for many years on research and development of such systems 
at Honeywell Inc. In a short note, Ohman gives a clear presentation of the 
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2.3 Sunnary 

The preliminary experiments carried out at this stage were critical to 
tne definition and continuation of the main research tasks. Not only was it 
true that the bulk of the work carried out later stemmed directly from this 
stage, but also most of tne instrumentation and measurement interpretation 
problems encountered when using infrared imaging systems in aerodynamic 
research, were crystalized during this stage. 
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there is a jump in the wall temperature, that can be explained by the fact 
that at the corner, the flow has no velocity components and therefore, the 
convective heat transfer there should be minimal. The second featjre was 
that at the point where the reattachment was believed to take place (* = 15 
cm), the temperature passed through a local minimum. It is a well known fact 
that the heat transfer coefficient in such a flow is the highest at the 
turbulent reattachment point (Merzkirch et al., 1988), thus explaining the 
minimum temperature value there. The third feature was that downstream from 
the reattachment point, the influence of the development of the turbulent 
boundary layer can be observed on the wall temperature distribution. The 
fourth and last feature of this experiment was that the assumed reattachment 
point occurred between six to seven step heights downstream of the step. For 
a Reynolds number based on the height of the step of approximately 17500, this 
result is in excellent agreement with the findings of Eaton and Johnston 
(1981). The results of this experiment showed that separated flows can be 
detected with infrared imaging systems by mapping the surface of interest. 
This conclusion led to the experiment described in Chap. 6 where the 
footprints of the airflow about a NACA 0012 airfoil were studied all the way 
from the fully attached regime at zero angle of attack up to the post stall 
regime. 

Before proceeding further with the analysis of the main series of 
experiments, it is essential to explain the operational principles of the 
infrared imaging system, in order to understand its limitations from an 
aerodynamic research point of view. 
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0.9x1. 2m ( 3 * x4 ' ) wind-tunnel. It is a low speed, return circuit wind tunnel 
driven by a 100 HP electric motor. The surface was warmed by a Constantan 
wire closely wrapped three times chordwise around the plate with a duct tape 
covering it in order to ensure uniform emittance and a smooth surface. The 
original idea was to produce a constant heat flux flat plate. However, it 
became very quickly apparent that the boundary condition produced is of 
constant power heating type. The difference between the two is that in the 
former case a constant heat rate is dissipated into the freestream through tne 
boundary layer, while in the latter the airflow, the surroundings and the 
substrate share the constant heat rate dissipated by the wires. The measured 
wall temperature distribution was considered physically realistic due to its 
similarity to the constant heat flux case (Gartenberg et al., 1987), except 
for the fact that after an initial raise with the square root of the distance 
from the leading edge it leveled off, instead of continuing the initial trend 
(Kays and Crowford, 1980, p. 151). It was concluded that this behavior was the 
result of an ever increasing thermal involvement of the substrate, as the 
laminar boundary layer thickens and the convective cooling becomes less 
effective. At this stage, it was felt that this experiment should be 
revisited with a wedge type leading edge that would permit comparison of the 
experimental with analytic results. Further details are given in Chap. 5. 

The last preliminary experiment was to monitor the wall temperature 
distribution for air flow over a rearward facing step to observe separated 
flow effects. The rearward facing step was produced by placing a right angle 
trapezoidal wedge on the flat plate previously mentioned, at 5 cm from the 
leading edge and parallel to it. The step height was 2 cm. The measured wall 
temperatures are reproduced in Fig. 2.3. Four features of the results were 
particularly encouraging. The first was that just at the foot of the step 
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of the mass per unit length with the specific heat) should be low in order to 
minimize the temperature time response to changes in local velocities. As 

already known, all these requirements led to the Chromel wire as the most 

suitable choice. On an absolute basis there are materials that fit better to 
the criteria established above. However, they fall into the category of rare 

metals, and thus their availability and price made their use in this case 
prohibitive. 

Another problem that arose during the heated wire experiments was the 

background radiation interference. Since (as we have seen previously) the 
wire is unresolved by the camera, photons originating from the background 

surface will hit the detector together with the photons from the target 
wire. The problem is aggravated with increasing background temperature and 
spatial nonuniformity. To prevent problems associated with background 
interference, its surface temperature should be kept at least uniform and 

preferrably low. All these factors indicate that the actual experiments 
should be run under the same background conditions at which the calibration 
was carried out. This is the reason why later, the laminar flow jet 

experiments (Chap. 4) were carried out against a uniform background made out 

of cardboard paper held at constant temperature by the air conditioning system 

of the laboratory. 

2.2 Surface Measurements for Boundary Layer Research 

A second series of preliminary experiments was concerned with detecting 
the laminar boundary layer development over a warmed flat plate. The 

experimental set-up is identical with that shown in Fig. 5.1 except for the 
fact that the leading edge was curved and not sharp as it should be in order 
to reproduce a Blasius type flow field. The experiments were carried out in a 
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Minimization of this factor is essential when velocity information is deduced 
from the temperature measurements. After tri al-and-error experimentations, it 
was decided that a 0.0762 mm (0.003") diameter Chromel wire would be the best 
compromise for the present task. At this point yet another problem was to be 
faced involving the wire surface emittance calibration. The Chromel wire has 
a highly polished surface with a typical emittance of 0.05. However, using 
this value to determine temperatures with the IR imaging system will result, 
as previously noted, in lower than actual values because the wire is 
unresolved by the camera. In order to get true temperature readings, an 
"apparent emittance" lower than the actual should be used. The calibration 
procedure is described in Chap. 4. In this context, it should be emphasized 
that for optically unresolved targets, the distance from the camera determines 
the fraction of the instantaneous field of view the target will cover on the 
focal plane. Therefore, emittance calibrations for unresolved targets are 
valid only for the specific distance at which they were carried out. 

The choice of the wire material should be made after a careful screening 
of its physical properties: the material should have a high electrical 

stivity in order to dissipate enough heat to raise its temperature and be 
visible with the IR imaging camera. Also, the temperature coefficient of 
electrical resistance should be low so that the wire will uniformly dissipate 
the heat lengthwise, even if the temperatures vary locally. The material 
should be chemically stable, otherwise corrosion may change the emittance, and 
thus accurate continuous temperature measurements may not be possible. The 
coefficient of thermal expansion should be low in order to prevent elongation 
due to the heating that may cause vibrations. The wire thermal conductivity 
should be low in order to prevent equalization of the temperature along the 
wire through conduction. Finally, the wire thermal capacitance (the product 
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Reynolds number. These heat transfer features are of course washed downstream 
with the flow causing the wake to be a zone of very intense turbulent mixing 
with peak values at the edges, tnus determining the temperature shape along 
the wire as shown in Fig. 2.2. 

An important factor that arose during this stage of the research was the 
cnoice of the wire material. From the flowfield po^nt of view, a thinner wire 
is more desirable. However, for lengths of wire herewith under consideration 
(tens of centimeters), a too thin wire may break easily. Also, a thin wire 
may not be resolved by toe camera i.e., its surface may not be large enough in 
order to be projected on the whole area of the camera detector. Therefore, 
the infrared imaging system wrongly assumes that the photons reaching tne 
detector come from an area which is equivalent at least to the instantaneous 
field of view of the camera for that focal distance, thus giving false lower 
temperature readings of the wire. (This aspect of the temperature 
measurements with the infrared imaging system will be addressed in greater 
detail in the next chapter). 

At this point, it seems that increasing the wire diameter may be a 
partial solution to the spatial resolution problem of the IR imaging camera 
but this is not completely true. A larger wire will disturb the flowfield, 
and also, as the wire diameter becomes larger, the spurious reflections from 
tne surroundings to the camera start to play an important role in the signal 
to noise ratio the camera receives. A separate set of experimental runs with 
a variety of wire sizes revealed that at 0.381 mm diameter (0.015") the signal 
to noise ratio becomes unacceptably large. For example, the presence of 
people in the laboratory could be detected through their body heat reflection 
from the wire into the camera. Later it was found that smaller diameters were 
also desirable in order to minimize the heat conduction along the wire. 
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experimental set-up was identical to that shown in Fig. 4.3 witn the 
addition that the wire was electrically heated. The measured temperature 
distribution along the wire was symmetrical with a minimum at the center of 
the jet, where the experimentally determined coefficient of heat transfer on 
the wire was in good agreeme/it with the value predicted by a standard 
correlation (Morgan, 1975). Although the result was considered physically 
correct, the opinion was that the heat-gun flow was not controllable enough in 
order to produce a flow field suitable for a thorough investigation. As a 
result, the conclusions drawn from this preliminary experiment led directly to 
the laminar flow jet experiment described in Chap. 4. 

In the second experiment, cylinder wakes were mapped usin^ a heated wire 
placed downstream, and again the temperature distribution was found with the 
infrared imaging system. The set-up as installed in the 0.9x1. 2 m (3 x4 ) low 
speed wind tunnel in which the experiments were carried out is shown in Fig. 
2.1. A 10 cm diameter cylinder was tested at a Reynolds number of 90,000, the 
wire being placed at 5 and 10 diameters downstream from the cylinder. A 
second cylinder, 2 cm in diameter, was tested at a Reynolds number of 18,000, 
the wire being placed at 8, 16, 24 and 32 diameters downstream. In all cases, 
the qualitative results were self-consistent and agreed with the physical 
expectations. A sample result is shown in Fig. 2.2 where the temperature is 
plotted as a function of the lengthwise wire coordinate. According to Schmidt 
and Werner (Schl ichting, 1968, p. 298), the heat transfer rate on the leeside 
of a cylinder, where the flow is separated, is equal to, or higher, than that 
on the windward side. The behaviour of the heat transfer rate in the 
separated flow area of the cylinder is dependent on the Reynolds number based 
on the cylinder diameter. In particular, just after the separation point, 
there is a jump in the heat transfer rate which is even more sensitive to the 
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Chapter 2 

THE PRELIMINARY EXPERIMENTS 

Since the main objective of this research was to extend the use of 
infrared imaging systems to a few important aerodynamic subjects of interest, 
the first step was to carry out initial experiments, where the suitability of 
the measurement technique could be easily judgeo according to tne results. 
The initial experiments and the results are described in detail ir. Gartenberg 
et al. (1987) and only the more interesting results are reproduced here. A 
review is made of the main ideas behind each experiment and of the way the 
results were evaluated. 


2.1 The Heated Wire Technique 

The first concept to be tested was the use of a thin, electrically heated 
wire as an air velocity detector. The idea was that by measuring temperatures 
along the wire while it is immersed in a flow with position dependent 
velocities, one may deduce the local velocities through Nusselt number 
correlations. This approach is useful when investigating shear flows, jets, 
boundary layer and confined flows, wakes etc. In such cases, the local heat 
transfer coefficient is mostly a function of the normal velocity component 
relative to the wire (Morgan, 1975); but even with this restriction, the 
method may be extremely useful to map flowfields unintrusi vely , since a 0.0762 
mm (0.003") slightly heated wire will not have any significant influence on 
the surrounding flow. Two experiments were performed in this category: 
mapping of an air jet and mapping of the wake behind a circular cylinder. 

In the first experiment, the 0.0762 mm (0.003") heated wire was placed at 
30 cm from the exit of a working heat-gun throwing cold air at 2.5 m/sec. The 
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The experimental work carried out under this program followed very 
closely the preceding stated purposes. Chapter 2 contains a brief review of 
an initial set of experiments that were aimed toward checking the feasibility 
of the objectives with regard to the heated wire concept and tracking the 
boundary layer behavior by surface thermography. Chapter 3 is a short 
description of the infrared imaging system, emphasizing the IR imaging camera 
operation principles. Chapter 4 describes the evaluation of the heated wire 
concept using a laminar circular jet as a test case. Chapter 5 presents the 
experimental results and an analysis of the laminar, thermal boundary layer, 
which develops over a flat plate, using surface thermography. Chapter 6 is a 
study of the boundary layer regimes over a NACA 0012 airfoil from zero angle 
of attack up to separation, employing also tufts that could be observed both 
visually and with the IR imaging system. Finally, in Chap. 7 the conclusion 
from this research are drawn and recommendations for further work are made. 
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With the realization that IR imaging systems can produce more useful data 
in aerodynamic research than has been generally appreciated, this work 
explores the extent to which these systems can be made useful in advancing the 
understanding of aerodynamic flow phenomena. As a result, the plan of 
research incorporated specific experimental objectives set out to produce an 
improved understanding regarding the detection of the low speed viscous 
interaction characteristics using IR imaging systems. Specifically, the 
following topics were identified to be pursued: 

!• Feasibility and validation studies of the heated wire concept. 

This method calls for placing a very thin and very long heated wire 
perpendicular to the velocity vector of an airflow of interest. By 
tracking the temperature distribution along the wire with the IR imaging 
camera, it is possible to capture the main features of the flowfield. 
Furthermore, by using appropriate Nusselt number correlations, the 
velocity distribution perpendicular to the wire may be deduced. 

2. Boundary layer development over a body of interest. 

Experiments and model analysis were planned so that by measuring with 
the infrared imaging system the wall temperatures of an actively heated 
body immersed in an airflow, the development of the thermal boundary layer 
can be tracked. 

3. Detection of boundary layer flow regimes over airfoils. 

The purpose in this experiment is to expand the applicability of 
infrared imaging systems to the detection of separation, beyond the 
detection of transition, e.g., Quast (1987). 
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For a constant Prandtl number, which is frequently the case in aerodynamics, 
Eq. (1.1) can be rewritten as 

£-2 - = Constant • * ■■-}" , 0<m<l (1«2) 

K a 

Let's assume a constant properties flows, as is the case in the subsonic 
regime where the temperature doesn't change appreciably throughout tne 

flowfield. Then, Eq. (1.2) can be expressed 

h = h (U, x) (1-3) 

where h is the local heat transfer coefficient, x is the coordinate location 
of the point of interest on the configuration surface, k a is the air 
conductivity, U is the local freestream velocity and v a is the air kinematic 
viscosity. This expression shows that the local heat transfer coefficient, h, 
between a surface and a surrounding flow is a function of the freestream flow 
velocity and the configuration geometry, x. if the air properties remain 
constant. On the other hand, the heat transfer coefficient h is defined by 

Newton's law: 

q" = h (T - T ) U* 4 ) 

'Iq W CO 

where q- Is the local heat flex, T„ is the local wall tercperatore and 
T is the freestream static temperature. 

ao 

Equations (1.3) and (1.4) are the basic formulae to be used for 
aerodynamic research using Infrared imaging systems. Measurement of the 
surface temperature distribution on the body of interest with the IR imaging 
system, and knowledge of the geometry, the freestream temperature and 
eventually the heat flux, may provide enough data to deduce useful information 
about the flow-body interaction. The thrust of the present research utilized 
this second approach. 
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concerned with the shuttle entry measurements concentrate on technical 
descriptions of the systems and the projected experiments hot give no 
resets. The other papers either go into lengthy descriptions of the infrared 
imaging sytem itself, or use them to produce artificially colored pictures, 
the quantitative data was either acquired or heavily backed-up by 
thermocouples. Each one of the references adds a little bit of information 
but none of them gives a complete picture, and this is the reason why every 
single reference that was found was listed. Of all the references, the most 

useful ones are those of Boy, an et al. (1978) from Arnold Air Force Station. 
Tennessee, and Quast (1987) from DFVLR in West Germany. 


1.4 The Present Work 

Considering the potential applications of infrared imaging systems in 

aerodynamic research, there are at least two distinct paths that may be 

followed. The first and the simplest one is to follow the Reynolds analogy 

between the heat transfer and the shin friction. Temperature measurements can 

performed in order to deduce Stanton numbers and subsequently skin friction 

coefficients. This approach works in fully developed laminar and turbulent 

flows, but it is not applicable in the region of the transition itself and in 
separated flows. 

The second approach is more fundamental and goes back to the basics of 
the convective heat transfer. As is well known, the dimensionless heat 
transfer coefficient, the Nusselt number (Nu), is usually expressed as a 
function of the Reynolds number (Re) and the Prandtl number (Pr): 


Nu = Constant • Re m Pr n 


( 1 . 1 ) 
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aircraft and some helicopters. The results show the very complicated airflow 
patterns that develop in ground effect and at low speeds. No follow-up of 
this experiment could be found. 


1.2 University Research 

Applicable references from universities around the world are very meager. 
Most of the investigations were done in the United States. It seems that 
Champagne et al. (1967) were first to use a microscope-pyrometer to make 
temperature measurements along a scaled-up hot-wire probe to determine the 
extent of the heat conduction to the supports. Meroney (1978) used an 
infrared imaging system to make heat transfer studies on model buildings in a 
wind tunnel. He proved the concept by estimating convection heat transfer 
coefficients but didn't proceed further because of lack of a suitable data 
processing system. Page et al. (1986) used an infrared imaging system to 
determine the temperature distribution in a heated jet at the stagnation point 
region. A similar research was done by Carlomagno and De Luca (1986) in Italy 
expanding the research to an array of jets and determining heat transfer 
coefficients from the temperature measurements. Finally, Spence (1986) used 
the infrared thermography to determine heat transfer coefficients on double 
wedge jet vanes used for missiles thrust vector control. 

1.3 Perspective on the Literature Survey 

The paucity of works documented in archival journals is a good indication 
regarding the immaturity of the infrared imaging method in aerodynamic 
research. Only three works, out of all those reviewed, have been documented 
in recognized scientific journals (e.g., Int. J. Heat Mass Transfer, J. Fluid 
Mechanics, AIAA J.). Furthermore, the relatively long reference list reviewed 
in this chapter is misleading with regard to its usefulness. The papers 
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that study in the light of the following experiment objectives: (a) determine 

the actual heating rate distribution on the lower surface, (b) the location of 
the boundary layer transition, and (c) location and extent of flow separation 
in front of the control surfaces. Green et al. (1983) give a detailed 
description of how the system works, together with a few representative 
results. However, no detailed presentation of the results or their processing 
followed up. Meanwhile, this research was cancelled, on financial grounds. 
In a second and completely separate experiment, the objective was set to be 
the measurement of surface temperatures on the leeside of the shuttle. In 
this case, an infrared imaging system is mounted in a pod on the tail of the 
Space Shuttle Columbia, from where it scans the orbiter's wings and fuselage. 
This experiment is of interest because of the highly vortical nature of the 
separated flow over the wings which would impinge on the leeside surfaces. 
The final objective is to compare the measurements with results of CFO codes 
(Anon, Aviation Week and Space Technology, Dec. 7, 1987). As in the previous 
case, no results have yet been published from this unique experiment. 

Contributing to the design of the pod, heat transfer measurements were 
performed by Nutt (1979). The radiation estimation calculations were done by 
Myrick and Kantsios (1982); the description of the infrared imaging system is 
given by Myrick and Throckmorton (1985); and, the overall experiment 

description is provided by Throckmorton et al. (1985). 

Two more noteworthy papers related to infrared imaging applications to 
flight testing are: one by Brandon, et al. (1988) from NASA Langley Research 

Center in which results relating to in-flight transition detection as well as 

mapping a vortex "footprint" on the airplane wing are described. Another 

paper is by Flaig (1977) from the Naval Air Systems Command, in which an 
infrared imaging system was used to map exhaust gases from a hovering VTOL 
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scanning direction is much poorer than across it. Tests were also carried out 
to determine the system spatial resolution using different size circular 
targets at the same temperature. Furthermore, an extensive measurement error 
analysis was carried out based on the work of Carter (1975 ). From all the 
aeronautical engineering groups using infrared imaging systems, the group at 
Arnold Air Force Station probably accomplished the most comprehensive analysis 
concerning the performance and limitations of this system in wind-tunnel 
testings. In two other papers, Stallings et al. (1979) and Stallings and 
Whetsel (1982) repeated some of the findings reported in the previous paper. 
Using the knowledge and facilities developed at Arnold Air Force Station, 
Martinez et al. (1978) performed heating experiments and measurements using an 
infrared imaging system on a 0.040 scale model of the Space Shuttle. ihese 
measurements served the Rockwell Corporation in the design of the thermal 
protection system of the orbiter and can be considered as one of the very few 
but nonetheless great contributions using an infrared imaging system in 
aeronautical engineering. Hender and Okabe (1983) made further use of that 
expertise- to carry out heating measurements on a wedge type model made of an 
elastomer bonded to an aluminium base. It was found that for the wedge 
configuration, the semi-infinite slab assumption breaks down when deducing 
heat transfer data. Although the outcome could be predicted a priori, its 
extent could hardly be so, and for such cases an infrared imaging system may 
be the cheapest and fastest means to obtain quantitative data. 

Two other applications relating to the Space Shuttle, concern mapping 
surface temperatures during the entry phase. In one of them, the temperatures 
on the windward side of the orbiters were measured with a telescopic infrared 
imaging system mounted on an airborne observatory. The feasibility analysis 
was performed by Swenson and Edsinger (1977). Later, Chocol (1979) reviewed 
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compared with the thermographies in order to determine the cooling 
effectiveness of the various configurations being tested (simple versus 
multirow holes). The research was performed in support of jet engines turbine 
design efforts, but no directly applicable results have been reported. 

1.1.5 U.S.A. 

Most of the infrared imaging systems applications in aeronautical 
engineering reported by national laboratories are linked to the Space Shuttle 
program. The beginning goes back to the initial design phases of the shuttle 
when the prediction of the heat transfer rates during the entry phase were of 
major concern. The first experiments were performed by Compton (1972) at NASA 
Ames. This experimental program was very similar to that of Thomann and Frisk 
(1967) in Sweden, except for the fact that extensive use of computers was made 

for both data aquisition and data reduction. However, the bulk of the 

experiments were carried out later on at Arnold Air Force Station in 
Tennessee. The initial study and set-up of the experimental system was done 
by Bynum et al. (1976). The study, which provided surface temperatures and 
heat transfer data, was still general in nature, the measurements being 
performed on a cone and a hemisphere at a Mach number of 8.0. Later, 

Stallings and Carver (1978) embarked on a program to make heat transfer 

measurements on Space Shuttle scale models, using an infrared imaging system 
and the phase-change paint method. Besides the obvious and original 
contribution to the Space Shuttle design, this group was the first to report 
the difficulties associated with infrared imaging system measurements. In a 
very extensive testing program, Boylan et al. (1978) evaluated the influence 
of the infrared imaging system optics on the overall performance. It was also 
showed that the system capability to map a step temperature change along the 
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wind tunnels that will allow transition studies at higher Reynolds numbers in 
ground test facilities. A first effort in this direction was to carry-out 
experiments at 240K. A short and succinct enumeration of the difficulties 
encountered when making IR measurements in cryogenic environments is given by 
Seraudie et al. (1988). A further contribution to the use of infrared imaging 
systems was made by Balageas and Ory (1985), who reviewed the previously 
mentioned semi -i nf i ni te slab assumption and showed in which cases it can lead 
to significant errors. To refine the calculation methods, a technique was 
developed that takes into account the curvature of the wing, its skin 
thickness, as well as the thermal boundary conditions at the back of the skin. 

1.1.3 West Germany 

According to a survey by Mordoff (1988), the Deutsche Forschungs-und 
Versuchsanstalt fur Luft-und Raumfahrt (DFVLR) conducts an extensive and 
systematic experimental program designed to prepare a data base of flow regime 
characteristics of wings and airfoils. The program includes both wind tunnel 
and flight tests (Quast, 1987), the infrared imaging system being the 
transition detection tool. Laminar flow data were recorded at Mach numbers 
varying up to 0.7, altitudes up to 10,000 m and sweep angles up to 23 
degrees. Further tests are planned to be carried out in order to expand the 
range of the parameters of interest. 

1.1.4 Japan 

A study concerning the film cooling effectiveness of injection from 
multirow holes was carried out at the National Aerospace Laboratory by Sasaki 
et al. (1979). In this study, an infrared imaging system was used as a 
visualization means for the spreading of the injection fluid. All the 
temperature measurements were done by thermocouples and their results were 
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results from the heat-conduction theory (e.g., Eckert and Drake, 1972, pp. 
168-176). This concept and the data reduction method was later adopted with 
little or no change at all by most of the groups working on this subject. 
However, Thomann and Frisk were not the first investigators to make the semi- 
infinite slab assmuption to deduce heat transfer rates from surface 
temperature measurements. This idea was borrowed from the phase change paint 
technique which later evolved into the thermographic phosphorescent paint 
technique. The interested reader can consult Jones and Hunt (1966) and 

Throckmorton (1972) regarding to the data reduction procedures originating 
from this technique. 

1.1.2 France 

The French at the Office National d’Etudes et de Recherches Aerospatiales 
(ONERA), were the first to realize the potential of infrared imaging systems 
to detect the laminar to turbulent boundary layer transition in the low 
subsonic regime (Bouchardy et al., 1983). This application is made possible 
by the fact that at the transition point there is an order of magnitude 
increase in the convective heat transfer rate. Thus, the transition line can 
be seen on an infrared imaging system display as a jump in the skin 
temperature going from low values in the laminar zone to high values in the 
turbulent zone, where the skin temperature is intially at the free stream 
temperature. Because at very low Mach numbers the temperature differences may 
not have been high enough to generate a good contrast on the IR thermography , 
the authors analyzed the raw data by applying digital image processing 
techniques. The effort was beneficial to other researchers as well, who made 
successful use of that capability (Schmitt and Chanetz, 1985). Lately, the 
effort at ONERA is concentrated on developing this method for use in cryogenic 
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flow, the boundary layer flow and the substrate thermal diffusion. This 
aspect was recognized a long time ago and led directly to the introduction of 
hot-wire and hot-film gages for boundary layer research. However, these 
sensors (like most of the others) can provide information that can be related 
only to their locations, and this is one of their shortcomings. It is the 
capability of infrared imaging systems to scan large areas in a dynamic 
fashion, and produce quantitative results related to the surface temperature 
that make them an attractive alternative aerodynamic research tool. 
Unfortunately, conservative thinking played a major role in the failure to 
recognize the full research potential of these systems. The following 
literature survey will point to the very few applications the ipfrared imaging 
systems have seen. The survey is divided in two parts: aeronautical research 
establishments and universities. The big gap in financial, manpower and time 
resources between the two had a clear impact on the span of the research 
output. 

1.1 Research at Aeronautical Establishments 


1.1.1 Sweden 

It seems that the first use of an infrared imaging system for aerodynamic 
research was made at the Aeronautical Research Institute of Sweden by Thomann 
and Frisk* (1967) which used an IR system to measure surface temperatures on a 
model in a hypersonic wind tunnel and to deduce heat transfer rates from the 
data. For reducing the data, the material beneath the model's skin was 
assumed to be a semi -infinite slab. This allowed to use well-known analytical 


*Complete citations are found in the alphabetical list of references section. 
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Chapter 1 
INTRODUCTION 

A quarter of a century after infrared imaging systems were introduced 
commercially, they are still rarely used in aerodynamic research. The 
perception of these systems as an exotic tool limited their application to a 
few specific tasks. This reality is rooted in the history and in the everyday 
practice of aeronautical engineering. Since the primal interest is on the 
forces acting on a given body, pressure measurements come out as the natural 
measurements choice. This traditional approach is very useful in design 
studies that are mainly concerned with forces and moments impressed on the 
body by the external flow. However, it has little to offer when it comes to 
understanding and analyzing the interaction between the body and the viscous 
fluid that takes place through the boundary layer. When viewed globally, the 
two views are, of course, interconnected because the geometry of the body and 
the freestream velocity will determine not only the pressure-dependent forces 
but also the character of the boundary layer whether laminar, turbulent, or 
separated. In the case of attached boundary layer flow, its regime, whether 
laminar or turbulent, will determine the magnitude of the skin friction 
drag. However, a separated boundary layer completely alters the external 
flowfield about the body and totally changes the pressure distribution on 
it. Thus, there is the need for a complimentary measurement technique that 
will offer insight into the whole picture of the fluid-body interaction. 

In the same way that the pressure distribution reflects the momentum 
exchange between a fluid and a body, the temperature distribution reflects the 
energy exchange between the two. The difference is that the energy exchange 
process is related to all of the interaction aspects related to the external 
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PREFACE 


This is a report of a research that was aimed at proving concepts, and as 
such was exploratory in nature. In order to keep the cost of the program 

within the assigned limits of the budget, some measurements were made with 
devices that were readily available in the laboratory. Some of these devices 
had an accuracy that may be judged as too coarse. At that time, the opinion 
was that as long as the measurement errors could be established quantita- 
tively, they sufficed for achieving the research purposes and paved the way 
for future investigators to apply the developed methods to their own research. 
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SUMMARY 


An experimental program was aimed at identifying areas in low speed 
aerodynamic research where infrared imaging systems can make significant 
contributions. Implementing a new technique, a long electrically heated wire 
was placed across a laminar jet. By measuring the temperature distribution 
along the wire with the IR imaging camera, the flow behavior was identified. 
Furthermore, using Nusselt number correlations , the velocity distribution 
could be deduced. The same approach was used to survey wakes behind cylinders 
in a wind-tunnel. This method is suited to investigate flows with position 
dependent velocities, e. g., boundary layers, confined flows, jets, wakes and 
shear layers. It was found that the IR imaging camera cannot accurately track 
high gradient temperature fields. A correction procedure was devised to 
account for this limitation. Other wind-tunnel experiments included tracking 
the development of the laminar boundary layer over a warmed flat plate by 
measuring the chordwise temperature distribution. This technique was applied 
also to the flow downstream *rom a rearward facing step. Finally, the IR 
imaging system was used to study boundary layer behavior over an airfoil at 
angles of attack from zero up to separation. The results were confirmed with 
tufts observable both visually and with the IR imaging camera. 
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maintained (or obtained) in order to compensate for the fact that the boundary 
layer thickens and heats up st teamwise. in this case, the constant heat 
generation rate is shared mainly by two participating media, the wooden 
substrate by conduction and the airflow by convection. Heat transfer by 

radiation is also taking part, but is insignificant. Since the air cooling by 

convection is more effective on the forward part of the plate, (due to the 
thin boundary layer), the heat transfer by conduction to the substrate is less 

than on the aft part of the plate, where the boundary layer is thick and the 

convective heat transfer is much poorer. The estimation of the boundary layer 
behavior will show that between the first and the last station of the first 

data group, it thickens by 901 (from 0.60 mm. to 1.15 mm) while with respect 

to the second data group, it continues to thicken but only by 30% (from 1.35 
mm. to 1.75 mm.) (see Appendix B, Eg. B7). 

Looking more attentively at Eig. 5.5, one may farther observe that the 

general trend of the slope of 0 vs. k 1 ' 2 is to decrease gradually from the 

leading edge downstream, at the same time that the overheat of the substrate 
is going up. This feature is an additional evidence to the ever increasing 
participation of the substrate due to the gradual decreasing effectiveness of 
the convective cooling, and it happens in spite of the fact that the wood 
qualifies as a heat insulator. 

The thermal energy analysis performed on the flat plate is summarized in 
Table 5.1. The convective heat transfer was calculated from the temperature 
regressions given in Fig. 5.5. It is a superposition of the heat flux as 
generated by a constant temperature distribution (Kays and Crawford, 1980, p. 

139) 

.1/3 Re 1/2 


Nu = 0.332 Pr 
x 


(5-2) 

: Pf r 

OF 'flip? 


Table 5.1 Mean values analysis of the heat dissipation modes on a flat plate heated at 
constant wall power generation placed in an airstream at zero angle of attack 
Re v = 1.233x10^ x (m). 
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(where Nu x is the local Nusselt number, Pr is the air Prandtl number and Re x 
is the local Reynolds number), and a constant heat flux as deduced from Eq. 
5.1. Thus, the regressions of Fig. 5.5 present two aspects of the convective 
heat transfer. The first is contiined in the point of intercept of the 
regression line with the 6-axis and is indicative of the level of overheat of 
the substrate that can oe accounted for through Eq. 5.2. The second is 
contained in the slope of the regression line, and is indicative of the 

variation in the chordwise temperature distribution, that in this case can be 
accounted through the constant heat flux model as given by Eq. 5.1. 

The quantitative evaluation of the experiment shows that the mean 
electrically generated power was 949 W/m 2 . Locally, along the centerline 
wire, where the temperature measurements were made, the value of the 

electrically generated power was 4184 W/m 2 (see Fig. 5.4). As the balance of 
energy analysis shows (see Table 5.1), the highest value of the mean heat flux 
convected into the air along that wire was 930 W/m 2 . This value is of the same 
magnitude to the solar heat flux at the earth surface. Therefore, this 

heating is thought to affect the boundary layer on the flat plate about the 
same as it occurs in reality over the wing of an airplane. 

5.4 The Infrared Imaging System as a CFD Codes Validation Tool 

The interaction between a solid body and a surrounding flowfield of 
different temperature has both a momentum and an energy exchange aspect. The 
former is responsible for the velocity distribution around the body while the 
latter is responsible for the heat, flux distribution over its surface. This 
heat flux will continuously change the temperature distribution of the body 
substrate and particularly of its surface. Thus, in general, even under 
steady flow with constant free-stream temperature conditions, the surface 

•. i\ ; 

(jf >*■> l* M 1 * 


temperature distribution of a body may remain time dependent for significant 
periods of time. 


Since the prediction of aerodynamic and aerothermal loads ay 
computational means is and will be one of the main research topics in 
aerodynamics, there is a clear need for an experimental method to provide data 
for comparison and validation of the CFD codes. Viewed from this perspective, 
the infrared imaging systems are ideal tools for such tasks since they do not 
require instrumented models, nor any specialized instrumentation conduits in 
the wind tunnel or on the tested aircraft. 

From a general stand point of view, the complete calculation of the 
interaction between an arbitrary body and the fluid flow requires the 
following steps: first, the Euler equation 



is solved for the external flow, yielding the flow velocity distribution 
around the body. Here, this equation is given in the incompressible flow 
formulation where, 1) is the air velocity vector, p is the air density, D/Dt is 
the total derivative operator, f is the body force and, p is the pressure. 

Next, the Navier-Stokes equation can be solved, perhaps in boundary layer 
approximation. This equation differs from the Euler equation by the inclusion 
of the viscous terms ^ V ti. Again, in the incompressible flow formulation 
it can be expressed as 

DU -* 2 -* 

P Of = F ' g r ad P + P a V U (5.4) 

2 

where 7 stands for the Laplacian operator and u a is the air viscosity. Now 
with the velocity distribution known around the body, one may solve the energy 
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equation for the temperature distribution in the flow field. Once again 
invoking the incompressiblity argument and assuming constant properties flow, 
thc energy equation as derive.l from its internal energy formulation can be 

expressed in terms of the air temperature field as 


r DT = R 
pc v Dt K < 


V 2 T + 


iT) 


'5.5) 


where c v is the air specific heat at constant volume, K a is the air tnermal 
conductivity, T is the air temperature and, * is the viscous dissipation 

function ( Schl i cht i ng , 1968, P- 254). 

ail the above mentioned equations with the addition of the continuity 
equation, can in principle be solved for the five parameters that umguely 
determine the flowfield behavior namely U, P and T. For the present purposes, 
the energy equation solution is of the main interest, since it .s the 
temperature gradient of the air in the boundary layer at the 
determines the heat flux into the body through the Fourier heat conduction law 


q' ' = - k a grad T 


(5.6) 


uhere „ the heat flux. In order to obtain the temperature distribution 

on the “body surface, it must be remembered that its surface is the solid 
boundary with respect to the fluid. Therefore, the heat conduction equation 


= a V 2 T 
dt 


[5.7) 


where a is the thermal diffusivity of the substrate, must he solved for the 
solid body in order to determine the surface temperature distribution. 

Boundary conditions do play a special role. For example, besides the 

f host <;nurces sinks and thermal radiation effects 
convective heat transfer, heat sources, 

may occur dt the solid boundary. 


The convection effect is more conveniently expressed through the Newton 
law 

= h (T w - U (5.8) 

rather than through the Fourier conduction law. Sometimes, the combination of 
the two laws is used to determine the convective heat transfer coefficient h, 
as a function of the difference between the solid wall temperature T w and the 
air freestream temperature T « 

CO 

The radiation heat transfer obeys the well known Stefan-Bo 1 tzmann 
equation 

= J T t (5.9) 

where a is the Stefan-Boltzmann constant and e is the radiation surface 
emittance factor. 

Schematically, this procedure is summarized in Fig. 5.7, which shows why, 
in general, this process is unsteady even for steady state flow conditions. 
The evolving nature of the temperature field of the substrate continuously 

changes the temperature gradient in the boundary layer and the resulting heat 

flux. 

Once the body surface temperature distribution is predicted, it is 

possible in principle to construct a synthetic temperature distribution of the 
body surface, as it is assumed that it will be produced by the infrared imag- 
ing system in a test run at identical conditions with those of the calcula- 

tion. Still to be remembered at this stage is the fact that the configuration 
geometry, the model material and its finish, and the IR imaging camera speci- 
fications, all these factors influence the final thermographic experimental 

result, and they should be taken into account when computing and synthesizing 
the predicted temperature distribution (Gartenberg and Roberts, 1988). 
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Hence, through a complete set of aerodynamic and heat transfer 

calculations , it is possible to come up with computed results mat are 
relatively easy to compare with experimental data obtained as thermographies . 
The whole concept is presented in Fig. 5.8. 


5.5 Summary 

The infrared imaging system is an effective, non-intrusi ve device for 
measurement of surface temperatures with position and time dependent 
temperatu re distributions. Its capability of mapping rapidly entire surfaces 
is unmatched by any other means existing today. Furthermore, for aerodynamic 
configurations, the surface temperatures provide important information about 
body interaction with the external airflow. This was basically shown for the 
case of a flat plate in the laminar boundary layer regime at constant wall 
heat generation rate. 

The concept of using infrared imaging systems in aerodynamic research as 
a direct CFD codes validation tool was introduced. In general, it can be 
stated that the surface temperatures of a body in a flowfield is time 
dependent even for steady state flows. 
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Chapter 6 


BOUNDARY LAYER REGIMES ON A NACA 0012 AIRFOIL 

The next step after the flat plate laminar boundary layer experiment was 
to expand the investigation to a two-dimensional NACA 0012 airfoil at various 
angles of attack from zero up to a post separation angle. Such an experiment 
offered the opportunity to validate the applicability of infrared imaging 
systems to aerodynamic research on a real aerodynamic geometry, and to provide 
some new information, for boundary layer flows extending all the way from 
fully attached laminar regime through transition to turbulen.ce, and up to 
separation. Furthermore, the extension of this method to three-dimensional 
geometries and higher Reynolds numbers flows can be easily implemented 
following the success of these experiments. 

6.1 The Experiment Viewed in Perspective 

The detection of boundary layer transition from laminar to turbulent 
regime or its separation from the surface of interest, is one of the most 

challenging subjects in aerodynamics. Almost any available experimental 
technique is used for this kind of studies: smoke visualization, laser 

Doppler velocimeters , tufts, chemical sublimation (acenaphtene) , oil flow, hot 
films and hot wires, phase change paints, liquid crystal paints, pressure 

measurements, optical techniques etc. With all these methods, the subject is 
far from being closed and the search for more rewarding techniques is still 

going on. As was already pointed out in the literature survey, infrared 

imaging systems have already been successfully used in the recent past to 
detect transition both in flight and in wind-tunnel testing. However, 

detection of separation with IR imaging systems is still an open subject 
because the heat transfer behavior itself in separated flows still resists 
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accurate prediction. Also, no vacation of tne » system finding 

was done in the very same experiments hy any independent means. Therefore. 
Crying to detect both the transition and the separation ,n the same experiment 
can enhance the usefulness of infrared imaging systems for aerodynamic 

research. 

6.2 The Experimental Set-Up 

As shown in Fig. 6.1, the experimental system consisted of a NACA 0012 
airfoil mounted vertically at one end of a rotating arm, at the other end of 
which was mounted the infrared imaging camera. The set-up was installed 
0. 9x1.2 m (3* x4‘ ) low speed, closed circuit wind tunnel. As the arm rotated 
around its center, the airfoil changed its angle of attack, but the camera 

remained in a fixed, perpendicular position relative to the airfoil chord. 

The angle of attack a. could be fixed within *0.5 degree of the desired 
value The airfoil had a chord of 0.3048 m (one foot), an aspect ratio of 1.5 
and it was clamped between two end plates measuring 35.5 x 7.62 cm <H"x3"> in 
order to reduce the three-dimensional flow effects of the tips. At zero angle 
of attack the airfoil was at a distance of about 10 cm from the wind-tunnel 
side wall, and the camera was similarly spaced near the opposite -all. While 
this arrangement did affect the airflow around the airfoil (by blockage and by 
wall proximity), the basic objective of the experiment, to detect the foot- 
prints of the boundary layer behavior with the IR imaging camera, was not 

hampered. The observations are reported as the, were made i.e., with no 

corrections to account for the wall proximity. In order to eliminate any 
possible spurious reflections from the airfoil to the camera, a 0.152 mm thick 
black matte paper was bonded to its surface with double-sided Scotch, tape 

. tho Airfoil with a 3 cm pitch. Figure 6.2 displays 
wrapped chordwise around the air 

the following features of the airfoil: 




yi 



Fig. 6.2 NACA 0012 airfoil with aluminum foil tufts 
in still air. 
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(1) Decimal subdivision of the upper surface arc length (further referred to 
as "s") was provided, where the lower and the upper divisions were marked 
respectively with highly reflective silver paint and ultrafine gold leaf 
to be visible on the IR imaging system display. A posteriori, it can oe 
concluded (see Fig. 6.4) that both the silver paint and the gold leaf 

markings performed equally well, and since the silver paint is easier to 

dpply it may b6 considorod as th6 pP 6 f 6 rabl 0 choice. 

(2) Above the arc subdivisions, four rows of aluminum foil tufts 0.0254mm 

thick and 2mm wide, observable both visually and with the infrared 

imaging system, were threaded through the paper in increments of one- 
tenth of the surface arc length (0.1s). 

(3) The surface was locally warmed by passing an electric current through a 
0.127mm (0.005 ) Constantan wire that had been wrapped chordwise around 
the airfoil underneath the black paper, at the following locations: 

(a) directly under the arc markings and under the tufts - to enhance 
their observability with the IR imaging camera by providing a "low 
temperature/low emittance" target against a "high temperature/high 
emittance background. 

(b) below the arc markings - to allow tracking the boundary layer 
development through chordwise temperature measurements on a smooth, 
uncontaminated airfoil surface. The idea to use this layout was 
suggested by the encouraging results obtained earlier with this 
method in observing the footprints of the development of a laminar 
boundary layer over a flat plate (Gartenberg et al., 1988). 
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The mean heating flax was approximately 500 «/■*. with no noticeable 
e tfect Oh the boundary layer based on the located and s,»ooth behuv.or of tne 
chordae surface temperature under the arc markings (it. M ■»««>■ The 
span-wise coordinates of these beating wires can be seen on fig. 5-3, marked 
by pencii on the black paper JU st below the arc markings. The uppermost and 

lowest lines demark the field of view of the camera. 

6.3 The Experiment Description 

The airfoil angle of attack was varied between 0" and 16". Photographs 
of the airfoil surface were taken under steady state conditions with a 36 - 
single lens reflex (SLR, camera, and of the >R imaging system display unit 

(oscilloscope face) using a Polaroid camera. In a separate se ^ 

i attark was varied between 0° and 16 in 

experiments, the airfoil angle of attack was 

4 . ,ic Anri the thermoqraphies were recorded 
increments of 2" at 2 seconds intervals, and the thermog 

on 5 i/4" diskettes by the system's dedicated microcomputer. 

The rationale behind the above steps was to validate by visual 
observations the tufts patterns taken with the IR imaging system, and to 
validate the temperature analysis performed on the microcomputer against 
observations of tbe airfoil surface and the tufts made with the IR imaging 
system display unit. In this way, future Investigations can rely entirely on 

observations performed with the IR imaging system. While it is true that e 

t t ran c it inn their contribution to the 
tufts themselves induce an ear y 

• coif PN/ident. Furthermore, transition induced by 
detection of separation is sel 

the tufts can be clearly seen on the IR imaging display. 

• v ,oro parried out at a nominal Reynolds number 
The wind tunnel experiments were c 

(Re c ) of 375,000, based on the airfoil chord. As the wind tunnel air heats up 
during operation, two opposite thermal effects take place concurrently on the 


airfoil surface: convective cooling of the electrically heated surfaced and 
convective heating of the cooler surfaces due to the heating of tne free 
stream wind-tunnel air above ambient conditions. The lateral conduction of 
heat was negligible, as it could be traced back on the thermograph i es . 

The heating-up of the wind-tunnel during its operation had a negative 
effect on the quantative interpretation of the computerized thermographies. 
In principle, the aim of this experiments is to look for surface temperature 
differences caused by changes in boundary layer regime. Qualitatively, tne 
change should be from lower surface temperatures under the laminar regime to 
higher temperatures under the turbulent (see next section). However, at trie 
same time there is an additional heating effect due to the heating of the 
freestream air in the wind tunnel. Thus, extracting temperature differences 
from thermographies taken at different angles of attack, loses its meaningful- 
ness due to the lack of a common freestream temperature. In other words, as 
the angle of attack of the airfoil is increased, there is an increase in 
surface temperatures due to changes in the boundary layer regime, and at the 
same time an additional change due to the air heat-up. Therefore, a 
posteriori, it can be remarked that the computerized thermographies could not 
be meaningfully interpreted and their use was discarded. At the same time, it 
is emphasized that the concept of analyzing transition and separation data 
through substraction of thermographies is sound and can work if the 
experiments are carried out in a thermally stable wind-tunnel. Also, this 
argument has nothing to do with regard to the determination of the transition 
line from the thermographies. On the contrary, the heating-up of the air only 
enhances the demarcation line of the transition. 


ORIGINAL FACE IS 
OF POOR QUALITY 
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6.4 Results and Discussion 

The results reported hereih are based on the observations made over 20 
different experiments. Figure 6.3 shows a sequence of the airfoil photographs 
taken of the airfoil surface with a 35 mm SLR camera, for angles of attack 
varying between zero to 13 degrees. These pictures partially document «e 
tufts behavior as they could be observed visually. The word "partially" is 
used here because continuous visual observation of the tufts offers more 
information than still photographs, especially with regard to the unsteady 
behavior of the flow. Figure 6.4 snows a similar sequence of pictures, this 
time as observed on the infrared imaging syst OT display. These thermographies 
were taken with the display unit set to the inverted mode which produces 
darker shades for higher temperatures. This choice was dictated by the fact 
that the tufts status can be better seen in the inverted mode than in tne 
normal mode. Thus, in Fig. 6.4, the darker shades represent higher apparent 
temperatures than the lighter shades. The pictures were taken with a 
Polaroid® camera and later reproduced with a 35mm SLR camera. These pictures 
prove that the tufts behavior is simultaneously observable visually and with 
the IR imaging camera, and therefore it is sufficient for future experiments 
to rely solely upon one set of photographies only, namely those taken from the 

IR imaging system display. 

Before going into the analysis of the main results, the facts on which 
the transition detection by thermography is based will be briefly reviewed. 
According to the Reynolds analogy 


c f 

St = o- or 


Mu 


Re Pr 2 


( 6 . 1 ) 


and under the assumption Pr = 1 
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35mm SLR camera photos. Air flows from left to right. 
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whore h is the local heat transfer coefficient Cf is the local skin friction 
coefficient, U is the freestream velocity, v d is the air kinematic viscosity 
and k d is its conductivity. According to Eq. (6.2), the local heat transfer 
coefficient is directly proportional to the local skin friction coefficient. 
Since during transition from laminar to turbulent regime, the skin friction 
coefficient jumps approximately by an order of magnitude, by Reynolds analogy 
the heat transfer coefficient also jumps by an order of magnitude. Thus, a 
clear jump in the wall temperature should be seen on the thermography as 
transition is induced. The temperature jump was shown to mark the transition 
region by prior workers (Bouchardy et al.; 1983, Quast, 1987), and can clearly 

be seen in Fig. 6.4, for example at a = 4° at the 0.35s location, on the arc 
markings . 

The main experimental results are discussed after giving a brief 
description of what is believed to happen on the NACA 0012 airfoil at relevant 
Reynolds numbers (375,000), as the angle of attack is increased from zero up 
to separation. This account is based on the visual tufts observation in these 
experiments as well as on the accounts of Carr et al. (1977, p. 3), McCullough 
and Gault (1975), and Thwaites (1987) pp. 200-205. At zero angle of attack, 
the flow is usually laminar over most of the airfoil. As the flow approaches 
the trailing edge, the laminar boundary layer separates with subsequent 
transition to turbulent flow. The expansion of the turbulent motion causes 
reattachment of the boundary layer as a turbulent one which remains attached 
up to the trailing edge. This localized region of laminar separated flow is 
called the laminar separation bubble. The air motion inside this bubble is 
circulatory with typical velocities of one order of magnitude lower than that 
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of the fre e stream velocity Rondel .od seller, 1988). The length of toe 
bubble is a few percent of the chord. As the angle of attack is Increased, 
the bubble and the transition region above it move upstream, thus increase 
she airfo.l area subjected to turbulent flow. At a certa.n point, tbe 
turbulent boundary layer fails to negate tbe adverse pressure gradients and 
part of the turbulent boundary layer separates, usualy near the tra 9 
edg e. This phenomenon is known as turbulent separation. As the angle of 
attack is further increased, the bubble moves upstream as does the point of 
turbulent separation, thus causing the gradient of the lift curve versus the 
angle of attack to decrease. As the bubble reaches the lead.ng edge, and the 
angle of attack is still increased, the separation point of the turbulent 
boundary layer continuously moves upstream, the bubble remain, ng relat.vely 
small and attached close to the leading edge. At a certain angle of attack 
(stall angle), the adverse pressure gradient can reach values where 
reattachement of the turbulent boundary layer ,s no longer possible, 
causes the bubble to burst over the entire surface of the airfoil from the 
leading to the trailing edge, causing what is known as "leading-edge bubble- 

bursting stal 1 

The main findings of the experiments (visual and IR imaging observations) 

, .. . Ku , a \ criission of some interesting points 

are summarized below, followed by 

worthy of elaboration: 

a = o Laminar flow up to 0.8s on the clean surface. 

Tufts induced transition at 0.6s 

a = 2 ° Transition at 0.7s on the clean airfoil (lower part), 

Tufts induced transition at 0.5s, 

High amplitude tufts fluttering after 0.7s 



" 4 ° Transition at 0.5s on the clean airfoil. 

Tufts induced transition at 0.3s, 

Higher amplitude tufts fluttering after 0.7s 

u _ 6 Transition at 0.2s on the clean airfoil. 

Tufts induced transition at 0.2s, 

Turbulent separation seen on tufts at 0.8s 

i = 8° Transition at 0.1s on all the airfoil, 

Turbulent separation at 0.7s 

a = 10° Fully developed turbulent flow all over the airfoil. 

Turbulent separation at 0.6s 

a = 12° Turbulent separation at 0.3s with flow buffeting, 

"Steady" separation at 0.5s 

a = 13° Bubble burst; complete separation. 

A first observation relates to the fact that the region of flow reversal 
due to separation could be very well seen in the tufts motion. Similarly, 
fully developed turbulent boundary layer flow could be identified in the tufts 
motion as a considerable increase in the flutter amplitude. As indicated by 
the tufts behavior, the turbulent boundary layer separation point was not 
stable with respect to time. Just after the separation had first been 
observed, oscillations of the separation line back and forth over at least 
one-tenth of the surface could be observed. Just prior to separation, these 
oscillations grew in amplitude up to one quarter of the surface. 

A second observation refers to the detection of the laminar separation 
bubble. In Fig. 6.4, this bubble can clearly be seen on the a = 6° thermo- 
graphy at 0.15s as a lighter shade curved strip. Since at the separation 
point the first derivative of the velocity, 5u/5y, is zero, so is the local 
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also the heat transfer coefficient h. 

sk^n friction coefficient c f , 5^ 

_ .n of lower somperdtures due to tne 

Therefore , the bubble appears as “■ 

l owe r local convective heating on a cool surface. Mow.ver, its observabi i J 
spends on a very fine tuning of the "then., -rtnoh of -display unit, 
wnich may he detrimental to the captore of other features, and may raise 

questions with regard to the reprohocibil . ty under any arbitrary conditions. 

r - ran rlparl v be seen in this 
Incidentally, the exteht of the bobble for « - u 

case to be about 10% of the upper surface. 

. firt that for steady state conditions 

A third observation relates to the 

we -in Fin 6 4 w^re taken, no substrate 
such as those under Which the photographs in F g. • 

, r ohi litv It seems that the heating 
heating was necessary for enhanced observability. 

of the wind tunnel, which in this case was KZK/.inute was more than enoug 

get a sharp image on the systems display. Also, the substrate chordw.se 

t-n have a well defined and localized effect (see ig 
heating can be seen to have a wen 

6 . 4 ). 

.. ct-aii anale of attack* f"Or the 
A fourth observation relates to the 

u hop _ tall a t 13° is somewhat premature. Hansman and 

Dresent test Reynolds number, stall 

* „ .tall angle of 14° for the NACA 0012 at a 
Craig (1987) for example report a st 9 

Reynolds number of 310,000. This premature stall could be expected considf i - 
in g the fact that the airfoil was placed near to and spanwise parallel t 

wind tunnel wall. This assembly geometry caused the a.rfoil to opera 

¥ effort nf increasing the actual angle of a a 
ground effect" which has a net effect of 

( Rae and Pope, 1984, p. «.)• ■« order to chect the correctness of s 

assumption, the airfoil was temporarily moved close to the center of th - 

tunnel . A few runs at the same Reynolds number showed repetitively that 
this position the airfoil did indeed staliat 14 degrees. 
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A closer look at Fig. 6.4 provides more detailed insight and helps 
correlate the observations with the flowfield features. For this purpose, the 

clean, undisturbed surface of the airfoil in the lower half of each 
thermography is noted. 


The thermographies for „ - 0, 2 and 4 degrees clearly show two heater 
zones. The first, starting at the leading edge is the effect of the conver- 
tive heating under the laminar boundary layer regime. The second, along the 
trailing edge and exhibiting a curved boundary shape is due to the effect of 
convective heating under the turbulent boundary layer regime. Just upstream 
of that curved boundary shape is the region where the transition from the 
laminar to turbulent regime takes place. As the angle of attack is increased 
from zero to two and from two to four degrees, the transition region advances 
upstream as expected. The next four thermographies from a = 6° to a = 12 ° 
display the flow regime where trailing edge turbulent separation is 
encountered. The main vortex resulting from the separation is located in that 
area, and since the turbulent mixing there is quite vigorous, the net effect 
is that heat is convected from the freesteam to the airfoil surface. In order 
to get a visual impression about the intensity of the turbulent mixing, the 
reader can consult Schl ichting (1968, p. 36) where a visualization of a 
separated flow about an airfoil is reproduced from Prandt 1-Tietjens. The 
value of the heat transfer coefficient in that area of separated flow is 
difficult to determine. However, measurements made by Schmidt and Wenner (see 
Schl ichting, 1968, p. 298) on cylinders in the supercritical range (Re = 
426,000) show a sixfold increase in the Nusselt number in the wake region 
immediately following the separation point, relative to the value immediately 
preceeding the separation point itself. This fact explains at least qualita- 
tively why there is a net heating effect in the region affected by the 
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trailing edge turbulent separated flow. The thermographic taken at 
and „ , 8' show the separation bubble moving upstream, and remaining small as 
long as the flow remains attached on the forward part of the airfoil (Carr 
et • a i., 1977 ). As the angle of attack is further increased to a * 10° 
and a = 12°, the stagnation point moves downward around the leading edge and 
followed by the bubble, thus getting out of the field of view of the iR 
imaging camera. As the angle of attack is increased to a= 13°, the stall 

angle is reached. However, the separation process is not immediate. It was 
observed on the tufts motion that at this angle of attack, the flow was 
extremely unsteady, the points of turbulent separation violently moving back 
and forth stream-wise before the flow separated itself completely from the 
airfoil. The dynamics of the separation appeared as if the separated flow 
tried to penetrate into the attached flow region and "peel it off" from the 
surface. The separation occurred usually spontaneously, without any change in 
the airfoil condition, after a very few seconds from the moment the airfoil 
reached a - 13°. The hystersis effect of the reattachment was also clearly 
seen, although no singular value for the reattachement angle could be fixed. 
It varied between 11 to 12 degrees. The reattachment itself was also preceded 
by extremely unsteady flow and strong buffeting. Usually, it could take the 
flow three to four failed attempts until it managed to reattach. These 
findings are documented in Fig. 6.5. Based on the tufts observation, it can 
be said only that the air motion in the separated flow regimes is that it was 

chaotically unsteady. 

The experimental results have been compared with predictions of the 
Eppler code (Eppler and Sommers, 1980). This computer program is accepted 
today as a good code with regard to inviscid and viscous flow calculations 
about airfoils. The code was run at a Reynolds number of 400,000. There is 
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Fig. 6.5 NACA0012 airfoil (a) at separation (b) at reattachement 

incidence. Photographs for each set taken five seconds apart 
under identical conditions. 
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good agreement between the code prediction and the experimental results 

regarding the stall angle. " 31 14 ^ 

experiment showed it to he 13 decrees near the wall and 14 degrees ,n the 

middle of the test section. At angles of attach helow the stall, the code 

consistent!, underpredicted the extent of the turbulent separation as well as 

its onset. The experiment detected the separation onset at six degrees (say 

seven degrees to account for the wall proximity) while the code predicted ,t 

at ten degrees angle of attach. The degree of agreement between the 

experimental results and the code prediction regarding the extent to which the 

airfoil surface is affected by the turbulent flow varies with the angle of 

attach, dp to four degrees, the code underpredicted the extent of turbulent 

fl „w approximately b, a factor of two (for partial support of this claim see 

Ladson 1983. p. 90, Fig. 30(a)). For six and eight degrees, the agreement 

starts to be very good (20% to 10% deviations) and from ten degrees and up the 

agreement is excellent. These results are summarized in Table 6.1. 

This chapter will be conciuded with a final remarh concerning the 
interpretation of the thermographies. Loohing at Fig. 6.4, one can observe 
that the temperature information contained in each individual thermography is 
not enough in order to allow a meaningful interpretation in itself. For 

example, the thermography tahen at 10 degrees angle of attack 
ambiguously interpreted as indicating laminar flow with transition 
turbulent flow at 0.76s, the actual case being transition to turbulent flow 
near the leading edge and separation at 0.76s (The tufts indicating the 
separation at 0.65s). However, a sequence of thermographies, covering the 

whole range of angles of attach, from zero up to separation, do provide enough 

data in order to evaluate the development of the boundary layer regime, from 

almost entirely laminar, through partially laminar, turbulent and separated 

and up to fully separated flow. 



Table 6.1 Transition an] separation studies over a NACA 0012 
airfoil at different angles of attack: 
experiment vs. Eppler code (Re c = 400,000) 


a 

, 1 
s turb 

thermography 

s turb 

Eppler 

Separated* 

tufts 

s separat.ed 

Eppler 

0° 

0.2 

0.4 

0.0 

0.0 

2° 

0.3 

0.6 

0.0 

0.0 

4° 

0.5 

0.8 

0.0 

0.0 

6° 

0.8 

1.0 

0.2 

0.0 

8° 

0.9 

1.0 

0.3 

0.0 

10° 

1.0 

1.0 

0.4 

0.1 

12° 

1.0 

1.0 

0.7 

0.5 

14° 

1.0 

1.0 

1.0 

0.9 


1 


In wall proximity. 
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6.5 Summary 

, he separation studies performed on the NACA 00:2 airfoil showed that 
aluminium foil tofts can be observed with infrared .using systems. This tact 
opens the possibility to use the infrared imaging system as a flow visualiza- 
tion and surface temperature measurement device at the same time. In 
agreement with other researchers, the laminar separation bubble and the 
transition were shown to be observable with the » imaging system. A rather 
intense and localized heat transfer process is seen to exist on the upper 
airfoil surface under the turbulent separated regime. Both the flow separa- 
tion and the reattachment were observed to be time-dependent phenomena. 
Regarding the complete separation dynamics, the observation of the tufts 
behavior indicated a "peeling off" effort of the separated flow region on that 
part of the boundary layer that is still attached. 


10 


Chapter 7 
CLOSURE 

7.1 Review of the Main Results 

The essence of the heated wire technique is to place a very thin and very 
long electrically heated wire in an airflow, the measurement of the 
temperature distribution along the wire with the IR imaging camera being a 
direct indication of the flow speed behavior. Furthermore, since the heat 
conduction and radiation are negligible compared to the convection, Nusselt 
number correlations can be used to calculate that airflow velocity component 
that is perpendicular to the wire. This experimental method may be considered 
as non-intrusi ve since the wire end supports, and the IR imaging camera, can 
be placed out of the flowfield. More importantly, the wire itself is very 
long, very thin and slightly heated, and therefore it doesn't really disturb 
the flowfield significantly. This method is not to be confused with the hot 
wire technique, where the support itself is immersed in the flow at the point 
of interest. 

The heated wire technique was evaluated using as a target flow a laminar 
flow circular jet. This type of flow produced large temperature gradients 
along the wire, up to 2200K/m. Starting at a temperature gradient of about 
350K/m, it was observed that the IR imaging camera has difficulties in 
producing accurate temperature results. Furthermore, the measurement error 
grew nonlinearly with the temperature gradient. It was determined that this 
phenomenon is associated with the very nature of this IR imaging camera and 
probably has to do with its minimum resolvable temperature difference 
performance (Ohman). The roots of this behavior could be tracked through the 
field of view and instantaneous field of view to the size and the 
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characteristics of the camera detector element itself. However, this actoal 
measurement error ,s always on the negative side of toe output, ,.e„ the 
readings are always less than the t-ue tempera, u-e. V,nce the measurement 

orror is well behaving with respect to the temperature gradient, it ts 

a nf | wp p-r-ar versus the temperature 

possible to produce a calibration curv 

gradients. Therefore, at least to a first order, the experimental data can be 
corrected to get a temperature distribution that is closer to the reality. 
This process can be repeated by trial and error, until the final corrected 
temperature distribution produces such a temperature gradient distribution 
that When subtracting the errors, corresponding to the adjusted temperature 
values, the original experimental distribution is recovered. 

The inverse approach of using the temperature readings and Nusselt number 
correlations to deduce the velocity distribution of the Jet was shown to worn 
in principle. With regard to this approach, it was shown how important U 
t0 work with accurate data, because of the great sensitivity of the velocity 
function to the temperature values. Unfortunately, error sources include not 
only the IR Imaging system, but also the heat transfer correlations 
themselves. As it was found out from the literature survey, this subject ,s 
far from being completely understood and therefore this problem should be 
looked upon from the broader aspect of the convective heat transfer 
discipline. The heated wire concept was also shown to work for more genera, 
flowf ields ..... cylinder wakes in w.nd-tunnel experiments, with very 

promising results. 

infrared imaging thermography has been previously used to determine 
convective heat transfer rates without paying attention to the development of 
boundary layer flow. This work has shown that it is possible to track the 
thermal boundary layer development over a surface by measuring the surface 
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temperature distribution with the JR imaging system. The method was proved on 
a flat plate at zero angle of attack, slightly heated chordwise at a constant 
power rate. In analyzing the experimental results, it is important in this 
case to account for all of the three modes of heat transfer: convection, 

conduction and radiation. Also, it is important to remember that due to the 
development of the thermal boundary layer over the surface of interest, the 
balance between the relative contributions of the three may change. m 
general, the surface temperature distribution was shown to be time dependent 
even for steady state flows, due to the evolving nature of the substrate 
thermal coupling with boundary layer heat exchange. Therefore, a sinyle 
thermography may not be enough in order to extract quantitative data about the 
flowfield, i.e., a time sequence is required. This method was also applied to 
the mapping of surface temperatures downstream a rearward facing step with 
very good results. The temperature profile correctly showed the lowest heat 
transfer rate to occur just after the step and the highest rate at the 
location of the turbulent boundary layer reattachment. Thereafter, the 

influence of the redevelopment of the turbulent boundary layer could be 
observed too. 


These experiments suggest that the infrared imaging system may serve as 
an experimental tool for validation of CFD codes since the predicted 
temperature distribution may be corrected for the directional emittance 
effects of the radiation. As a result, the calculated surface temperature 
distribution as assumed to be produced by the IR imaging camera may be 
directly compared with the results of an actual thermography resulting from an 
experiment run under identical conditions. Generally, it is easier to run 
first the experiment and then to run the codes, since a precise tuning of the 
experimental condition nay be sometimes an impossible task. 
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t )f detect i no lew .need boundary layer transition over 
The concept ot oetecun. 

- . K M,t fi'/ti years old. The present 
airfoils using IR imaging systems b 

f ^n*il ir ion to turbulent repartition and laminar 
research extended the range of -Pi 

„ n u was show., to bo possible to detect the 
complete separation as well. 

v ari„us boundary layer regimes, I,.. 1-1"-. " 

and separated flow, all on one thermography imaging a MA 0012 amtoi, at 

of atuck j„ ground effect at a- chord Reynolds number o 
degrees angle of attao. 3 

375000. The development of tufts that can be observed both visually an w, 
the IR imaging camera greatly enhances the capability of this too, in 
aerodynamic research, because the method may serve simultaneously in t e 

temperature measurment and the visualization roles. Rased on the observation 

that the flow separation and reattachment 
of the tufts behavior, it was seen that the now s p 

are phenomena that develop in time. The flow prior to both of them ,s high y 
unstable and unsteady, a phenomenon commonly known as buffeting. Kith regar 
to the mechanism of tne total flow separation, also known as leading edge 
laminar bubble burst. It was observed that there is a pulsating behavior o 

the turbulent separated boundary layer flow back and forth chordwise, 

• • n nnrtion of the attached flow from the 
tries to "peel-off" the remaining portion or 

airfoil surface. 

finally, the reader is reminded that the infrared imaging camera is 

basically a device sensitive to the number of photons originating from the 

instantaneous field of view which reach the detector. H is through this 

mechanism that the IR imaging camera output is sensitive to the sur ace 

f Therefore in extracting quantitative data from 

temperature. Thererore, 

. nt to .-emember the factors that may influence the 
thermographies it is important 

i and reflections, directional emittance 

readings, namely emittance values and ret lectio 

effects, depth of field and proper focusing etc. 


7*2 Recommendations for Future Work 


Since 

nature, it 
should be 
rather than 


■nest of fie work performed during this program was new 
is f e , t that each one of the various experiments reported 
revisited, concentrating on one at a time, stressing the 
the span of the research. 


in i i b 
here \ n 
depth 


The heated wire technique can be applied to virtually any kind of 
subsonic flow of interest. Boundary layers (laminar or turbulent), confined 
flows, free shear layers, wakes, jets etc. can be studied using this method. 
Of course, the first item to consider is the temperature gradient induced 
error. The precise behavior of this function is required before starting any 
new experiment. It should also be remembered that each calibration is valid 
only for the specific conditions under which it was performed. One question 
to be answered is whether there is a possibility to produce one universal 
curve that will show the measurement error as a function of the temperature 
change per instantaneous field of view rather than per unit length at the 
target. For a resolved target, this task may be easy to accomplish, and with 
some effort it may be possible even to determine a correction factor that will 
account for various degrees of unresolved targets. One method to check this 
possibility is to map the same target with various IR imaging camera objec- 
tives, thus changing both the field of view and its ratio to the instantaneous 
field of view. Referring to the heated wire experiments, it may be useful to 
know that the infrared imaging camera may be transformed to a one-dimensional 
scanning device by stopping the rotation of one of its scanning prisms. By 
stopping the horizontal prism for example, it is possible to get a scanning 
rate of 2500 Hz of a certain line, and this raises the possibility of making 
frequency response experiments using the heated wire technique. 
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• c .,rnM -ns with the inverse approach of deducing 
One of the more serious prob Uis w 

no. velocities fro* temperature measurements is that the heat transfer rate 
„ greatly influence* hy the free stream turbulence level ^chUchting, 1*8. 

299). Thus, it is possible to use a Nusselt number correlation base* on 
the wire diameter an* deduce results based on the assumption that the flo 

laminar, while in reality the flow is turbulent and henaving differently wit 

, u„ Ine of the problems with the laminar flow 
respect to the heat transfer mode. One of P 

get experiment was that at a diameter based pipe Reynolds number of abou 

100 0, the flow starts its transition to turbulence, bub fully accomplishes , 

only at a Reynolds number of about 3000 (Blevin, 1984, p. 231). There ore, 

change in the flow velocity, could produce different results not only because 

inwpi nf the turbulence was 

the velocity was changed, but also because 
changed. 

In general, it is felt that this method is very suitable for heat 
transfer research, such as the influence of the freestream turbulence leve , 
or to study the process at low Reynolds or Grashof numbers. 

With regard to using surface thermograph, for boundary layer research, it 
is very important to have agreement between the expected surface thermal 
response and the phenomenon investigated. A variety of possibilities are 

heated prior to, or during the whole duration of the experiment ,n order o 

follow the development of the thermal boundary layer. 

lofort-iin experiments should really concentrate 
The boundary layer regime detection exper 

on the separation signature. The infrared "visible" tuns method should be 
qualified also on flight tests, although some trouble can be expected ,n e 
air because the tufts wil, be unresolved by the camera. This fact, combined 
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with their low emittance may cause the tufts to become "invisible" to the IR 
imaging camera during flight tests. Regarding the flight tests, it should be 
emphasized that the ultimate goal in developing this experimental tool is that 
by scanning an airplane wing in flight, the test engineer will be able to 
determine the nature of boundary layer behavior at a specific station i.e., 
laminar, turbulent or separated. 


7.3 Conclusions 

The purpose of this study was to explore new applications in which 
infrared imaging systems could contribute to and enhance aerodynamic 
research. Two different experimental approaches have been studied: the 

heated wire technique and the surface thermography. The heated wire technique 
is a completely new approach to the study of flows with position dependent 
velocities such as jets, boundary layers, shear flows, wakes etc. A 
significant finding enabled the measurement error to be represented as a 
function of the temperature gradient, and provided a method of correction for 
the experimental results in order to get the actual temperature distribution. 
The use of surface thermography was extended to research on boundary layer 
flow detection spanning all the various flow regimes; laminar, transitional , 
turbulent and separated, and for the study of separated flows following a 
backward facing step. In order to confirm the thermographic results, a 

visualization method was developed using tufts detectable by the IR imaging 
camera, thus making this system a simultaneous, temperature measurement as 
well as a flow visualization tool. 
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Appendix A 

APPENDIX TO CHAPTER 4: THE LAMINAR FLOU JET EXPERIMENT 

Ai: Error Calculations 

These calculations ignore the errors that, may originate in the 
difficulty to track nigh temperature gradients. 

1. Wire temperature read-outs from Disco 

(a) Due to uncertainty of the ambient temperature 

A j = 0.05°C AT w = 0.9°C 
d 

1 - Amiracv (through the transmittance factor in 

(b) Due to emittance input accuracy \ y 

Disco II). The transmittance last significant digit is !%• A 
transmittance difference of 0.5% implies a 
at = 0.2°C 

W 

(C) Due to temperature induced changes in the surface emittance - unknown 
in nature 

(d) Due to improper focusings - practically zero. 

(e) Due to the AGA infrared camera accuracy &T W - 0.2 C 


AT 


= ’^oTg 1 2 +0.2^ + 0.2 = 0.9°C 

w 


Therefore: 


2. Moan jet velocity assumed at a certain point (pixel) 


id? 

(a) Due to uncertainty in total pressure and temperature measurements 
ahead of the sonic nozzle: 

Au = 0.09 m/sec. 

Note: The total temperature induced error is 2 orders of magnitudes 
lower than that induced by the pressure error. The pressure error 
has three contributions: 

i) reading accuracy + 1 psig 

ii) pressure fluctuations + 1 psig 

iii) Gage accuracy assumed a 3% i.e., + 0.9 psig. 
Altogether this accounts for + 1.68 psig. 

(b) Due to uncertainty in the nozzle's discharge coefficient 1%: 

Au = 0.02 m/sec 

(c) Due to the camera's resolution (translates into half a pixel on the 
computer's screen) 

Au = 0.36 m/ sec 
Therefore : 

Au = ^(0.09)^ 7 (0.02) 2 + (0.36) 2 = 0.37 m/sec 

3. (a) Error in the calculated temperature profile due to velocity 

uncertainities : ±1 % from the absolute temperature in deg K. 

(b) Error in deduced velocity from temperature measurements: ±15% from 
the local velocity. 
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The error in the calculated temperature was estimated with data from 
experiment #3 with Au ~ 9.37 m/sec. The error in the deduced velocity 
was estimated in the same way with aT - 0.9K. 

4. The relatively high errors encountered in this experiment should be 

credited in their greatest part to the purpose of the experiment, that was 
to evaluate to what extent the IR imaging system can detect relatively 
nigh temperature changes that occur over a small portion of the camera 
field of view. When viewed in this light, the system's performance is 

quite satisfactory. This fact can be better understood by reconsidering 
the magnitude of the error induced by the uncertainty in determining the 
exact location of the measured spot (the "half pixel" uncertainty, see 
item above). If, for example, the jet would have been mapped along the 
whole width of the field of view of the IR imaging camera, this error 

could be reduced by an order of magnitude. This of course, would require 
a different objective lens. 

5. The error in determining the velocity distribution in the core flow 

of the jet varied between 2.8% on the centerline and 12.5% near the 
wall. According to the manufeturer ' s specifications, the Pitot static 
tube may induce an error of +2%. The rest of the error originates in the 

accuracy of reading of the difference between the water columns in the U- 

shaped manometer (+0.5 mm). 

A2: Miscellaneous Calculations 

(1) Radiation losses estimation: 


a = e = 0.05 



:.’4 

E = en(T 4 - T 4 )A = 0.05x5.668xi0” 8 (373 4 -29J 4 )* n ♦ 7.62xlO' 5 xl .0 

W d 

= 0.008 w/m lenyth of wire. 

(2) Conduction losses estimation (per pixel distance) for 
cjelec = 9.33 w/m (maximum heating rate): 

1 pixel at 0.5 m focus = 1.154 mm. 

<3 = - kA[grad T - grad T ] • — - - — 

a ^ ^ ^ uA i 2 

T(°C) = 78.1 - 10.2 eos[ -r~r—~~r w] ; ( r 2 = 98.3%; a = IK) 

1 Z . 6^4 


grad T = 10.2 (- 


12.694x10 


— t) S1 n ^ To~T57 11 -1 ( K / ,n ) 
-3 "12.694 J 


AQ = 50 • - [7 .62xl0 _5 ] xl0.2 

4 -3 

conduction 12.694x10 

max 


[0.989-0.910 


= 3.97x10” w/m length of wine 


(3) The free convection estimations were based on Kays and Crawford (1980) p. 
328: 

Nu » 0.518 [1 - (^|22)°* 6 f S/12 (Gr 0 Pr) 1 ^ (A.l) 

with the correction for thick boundary layer compared to the cylinder 
diameter 


Nu 


2 

In (1 + 2/Nu) 


(A. 2) 
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The transcendental Eqs. for AT assuming T^. = 323K is 


at 


887 


2 

\ 


1 n 


1 + 


2 

, ' 0.25 

0.029(aT) 


t A-3) 


where i the electrical current in Amperes. 


(4) The coefficient of determination r 2 (adjusted for degrees of freedom) is 
given by: 

? . Sum of squares of t he residuals /(n-p) (A. 4) 

r = 1 " Sum of squares of the total values /( n-1 ) 

where n is the degrees of freedom and p is the number of the coefficients 

fit in the regression equation. 


( 5 ) Calculation of the predicted wire temperatures in Figs. 4.6 and 4.7 and 
of the deduced velocity profiles in Figs. 4.10, 4.11 and 4.12. 

The correlation (Eq. 4.1) 


Nu = 0.795 Re 


0.384 


(A. 5) 


can be explicitly written 


{Li = 0.795 (— ) 

V 


0.384 


(A. 6) 


where h is the heat transfer coefficient, d is the wire diameter, k a is 

the conductivity of the air, U is the local velocity of the air 

and v is the kinematic viscosity of the air. 
a 


Defining AT = T wi 


re 


T air 


h “ XvT 


and 


where q is the heat 
the ci rcumferential 

at = • 


generated per unit 

area of the wire per 

• v 0.384 

3 r_il u 

0 . 7 95 ti k d ' IJ 


a 


length 

unit 

-0.384 


of wire (q = i ? R) and A is 
length (A = id), one gets 

(A. 7) 


This formula was used as is to deduce air velocities from the wire 
temperature measurements (Figs. 4.10-4.12). To predict tne wire 

temperatures in Figs. 4.6 and 4.7, the velocity profile was assumed to 
have a parabolic distribution: 

U < r) ■ U max h-(f) 2 J (A.S) 

where U max is the centerline velocity of the air in the pipe and has 
twice the value of the mean velocity. 

A3: Approximating Functions for the Wire 

Temperature Distributions 

Assuming a cosine function approximation for the temperature variation 
along the wire, with the maximum temperatures occurring at X = 0,4 (0 and 2 n 
radians respectively), and the minimum at X = 4/2 ( n/2 radians), the tempera- 

ture profile can be approximated by: 

T(x) = 0 . 5 T ( T + T - ) + (T - T - )cos ^2-1 (A 91 

' max min' v max min' 4 J 

In Figs. 4.6 and 4.7, zero radians correspond to pixel #-11, n radians to 
pixel #0, and 2n radians to pixel #11. Therefore 

grad T(x) = -(T - T . ) -2 sin -^2- (A. 10) 

max min' 4 4 1 


The evaluation for the four experiments is given in Table A.l below. 
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Table A.l 


Approximating functions for the wire 
distribution and temperature gradient 


emperature 




H = 0.0254 m 

Ax = 0.001154 m (1 pixel) 


Appendix R 
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APPENDIX TO CHAPTER 5: THE FIAT PLATE EXPERIMENT 

Bl: Error Calculations 

1.1 Surface temperature read-outs from Oisco 11 (software package for 1R rtdU 
reduction) 

(a) Due to uncertainty of the ambient temperature, practically 

zero (AT << aT ). 

® w 

(b) Due to uncertainty of emittance value. 

Ae = 0.02 AT - 0.6°C 

(c) Due to temperature induced changes in the surface emittance - 
unknown in nature. 

(d) Due to unproper focusing - practically zero. 

(e) Due to the AGA camera resolution 

AT ( = 0.2°C 
w 

(f) Due to error in location 

AT = 0.1 5°C 
w 

AT W = + 0.2^ + 0.15 2 = 0.65°C 


Therefore, 
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1.2 Error of the slope of •> vs x 1//2 (mean value,; 

d:-) 2x 0 .65 


' s “ 1/2, 0.4 - 0.2 

d x 


°C 

6.5 — 
m 


1.3 Error of velocity measurement 
Estimated at 2% 


All = 21 < 0.02 = 0.42 m/sec 


1.4 Error in the (calculated) convected heat flux (worst case); 


1/3 ,1/2 
de k Pr U 

T/27 1/2 


0 d(x ) 2.2019 v 


A q„ = 


t Pr 1/3 {u l/2 . 4 


- A'J } 


2.2019 v 


1/2 


0.02709 • 0.709 


1/3 


d (x 1/2 ) y d(x 1/2 ) 2 U 1/2 


1/2 50.07 

f 20. 9 • 6.5 + rrr * °- 42 ) 


2.2019 • (16.95 x 10 


"" 2 -20.9 1/2 


= 2.66 (29.72 + 2.3) = 85.2 ^ 

m 

Thi.s is about one fifth of the mean convected heat-flux value and it 
comes almost entirely from the error in temperature measurements, which 
originate in the uncertainty of the emittance value (see items 1.1 and 1.2 
above). This is a well known difficulty and is discussed in Sparrow and Cess 

(1978). 

B2: Miscellaneous Calculations 


1 . 


Determination of the temperature distribution along a flat plate in the 
laminar boundary layer regime with constant heat flux, equation (5.1). 


Starting from equation (9-39) in Kays and Crawford (1980) 


L3U 


T (x) - T = Pr" 1/3 Re' 1/2 f X ri - (V /4 r 2/3 q‘ ‘(p) d- (B.l) 

o « k -o k x J o 

• r 1 
for 1 ( r, ) = const, and with y = — ; dy = — d£ the integral becomes 


r 1 3/4 r 2/3 , 

x f 1 1 - y ) dy 

o 


(B.2) 


3/4 3 - 1/4 4 • 1/4 

Now z = y and dz = — y dy -*• dy = — dz (y ) and the integral 

H J 

becomes 


f l f . 2/3 4 1/3 

x J (1 - z) • - z dz 

o ' 3 


, . ,1 m-1 . , . n-1 , , r(m) • r(n) 

Since p (m,n) = / z (1-z) dz = p (m,n) = - 

1 o 1 r(m+n) 

the integral 


4x ,1 . i -1 1 \ -1 1 

— / (1 - z) 3 z 3 dz n = 

3 J o v ; 


3 m = 1 3 


(8.3) 
we get for 

( B.4 ) 


Since r(n + 1) = n r(n) f( j) = 3r (l y) , we get 


4x ^ ¥ r(1 T 1 

3 r (' f) 


4x 3 [hi j)f 
3 ,, 2x 

r(i 3) 


4x 3(0.89338)' 
3 0.90330 


3.5343X 


( B.5 ) 


Inserting this value for the integral in (B.l) and defining 
G(x) = T (x) - one gets 


0(x) = 2.2019 ^ x Pr" 1/3 Re" 1/2 
k 


(B.6) 


which is precisely equation (5.1). 


2. Calculation of the velocity boundary layer thickness. 


Following Schlichting, (1968) p. 26, the laminar velocity boundary 
layer thickness can be well approximated by 



131 


- 1/2 

6 = 5 Re x l/ x or 

, v , 1/2 1/2 

6 - 5H x 


The numerical value for this experiment was calculated accord.no to 


1/2 


6(m) = 0.004375 [x(m)] 


(B.7 ) 
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